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Chapter 1 

General introduction 



1.1 Preface 

This thesis deals with the study of the structure and the interactions of the fundamental constituents 
of matter. We arrived at the end of the twentieth century describing the known fundamental prop- 
erties of nature in terms of quantum field theories (for the electromagnetic, weak and strong nuclear 
forces) and general relativity (for the gravitational force). The Standard Model (SM) of the fun- 
damental interactions of nature comprises quantum field theories for the electromagnetic (quantum 
electrodynamics, QED) and weak interactions (which are unified in the so-called electro-weak the- 
ory) and for the strong interactions (quantum chromodynamics, QCD). This SM is supplemented 
by the classical (not-quantum) theory of gravitation (general relativity). All the experiments that 
has been performed in accelerators (to study the basic constituents of nature), up to now, are con- 
sistent with this framework. This has led us to the beginning of the twenty-first century waiting for 
the next generation of experiments to unleash physics beyond that well established theories. Great 
expectations have been put on the Large Hadron Collider (LHC) actually being built at CERN 
and scheduled to start being operative in 2007. It is hoped that LHC will open us the way to 
new physics. Either by discovering the Higgs particle (the last yet-unseen particle in the SM) and 
triggering this way the discovery of new particles beyond the SM or by showing that there is not 
Higgs particle at all, which would demand a whole new theoretical framework for the explanation 
of the fundamental interactions in nature 1 . Possible extensions of the SM have been widely studied 
during the last years. The expectation is that those effects will show up in that new generation of 
experiments. Obviously accelerator Earth-based experiments are not the only source of information 
for new physics, looking at the sky and at the information that comes from it (highly energetic 
particles, cosmic backgrounds...) is also a widely studied and great option. 

But in this way of finding the next up-to-now-most-fundamcntal theory of nature we do not want 
to lose the ability to use it to make concise predictions for as many processes as possible and we also 
want to be able to understand how the previous theory can be obtained from it, in an unambiguous 
way. It is obviously the dream of all physicists to obtain an unified framework for explaining the four 
known interactions of nature. But not at the price of having a theory that can explain everything 
but is so complicated that does not explain anything. In that sense, constructing a new theory is as 
important as developing appropriate tools for its use. As mentioned, this is true in a two-fold way, 
we should be able to understand how the previous theory can be derived from the new one and we 

1 Let us do not worry much about the scaring possibility that LHC finds the Higgs, closes the SM, and shows that 
there are no new physics effects at any scale we arc capable to reach in accelerators. Although this is possible, it is, 
of course, extremely undesirable. 
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also should be able to precisely predict as many observables as possible (to see if the observations 
can really be accommodated in our theory or if new effects are needed). Let us end this preface to 
the thesis with a little joke. According to what have been said here, the title of the seminar that 
would suppose the end of (theoretical) physics is not: M-theory: a unification of all the interactions 
in nature, but rather: How to obtain the metabolism of a cow from M-theory. 



In the thesis we will focus in the study of systems involving the strong interacting sector of the 
Standard Model (SM). The piece of the SM which describes the strong interactions is known as 
Quantum ChromoDynamics (QCD). QCD is a non-abelian ST/ (3) quantum field theory, that de- 
scribes the interactions between quarks and gluons. Its Lagrangian is extremely simple and it is 
given by 



In that equation qi are the quark fields, igG^ v = [D^, D u ], with 0^=8^+ igA^, are the gluon 
fields and Nf is the total number of quark flavors. QCD enjoys the properties of confinement and 
asymptotic freedom. The strong coupling constant becomes large at small energies and tends to 
zero at large energies. At large energies quarks and gluons behave as free particles. Whereas at low 
energies they are confined inside color singlet hadrons. QCD develops an intrinsic scale Aqcd at 
low energies, which gives the main contribution to the mass of most hadrons. Aqcd can be thought 
of in several, slightly different, ways, but it is basically the scale where the strong coupling constant 
becomes order one (and perturbative calculations in a s are no longer reliable). It can be thought as 
some scale around the mass of the proton. The presence of this intrinsic scale Aqc d and the related 
fact that the spectrum of the theory consists of color singlet hadronic states, causes that direct QCD 
calculations may be very complicated (if not impossible) for many physical systems of interest. The 
techniques known as Effective Field Theories (EFT) will help us in this task. 

In general in quantum field theory, the study of any process which involve more than one relevant 
physical scale is complicated. The calculations (and integrals) that will appear can become very 
cumbersome if more than one scale enters in them. The idea will be then to construct a new theory 
(the effective theory) derived from the fundamental one, in such a way that it just involves the 
relevant degrees of freedom for the particular energy regime we are interested in. The general idea 
underlying the EFT techniques is simply the following one: to describe physics in a particular energy 
region we do not need to know the detailed dynamics of the other regions. Obviously this is a very 
well known and commonly believed fact. For instance, to describe a chemical reaction one does 
not need to know about the quantum electrodynamical interaction between the photons and the 
electrons, but rather a model of the atom with a nucleus and orbiting electrons is more adequate. 
And one does not need to use this atomic model to describe a macroscopic biological process. The 
implementation of this, commonly known, idea in the framework of quantum field theories is what 
is known under the generic name of Effective Field Theories. As mentioned before, those techniques 
are specially useful for processes involving the strong interacting sector of the SM, which is in what 
this thesis focus. The process of constructing an EFT comprises the following general steps. First 
one has to identify the relevant degrees of freedom for the process one is interested in. Then one 
should make use of the symmetries that are present for the problem at hand and finally any hierarchy 
of energy scales should be exploited. It is important to notice that the EFT is constructed in such a 
way that it gives equivalent physical results (equivalent to the fundamental theory) in its region of 
validity. We are not constructing a model for the process we want to study, but rigorously deriving 
the desired results from the fundamental theory, in a well controlled expansion. 



1.2 Effective Field Theories 




(i.i) 
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More concretely, in this thesis wc will focus in the study of systems involving heavy quarks. As it 
is well known, there are six flavors of quarks in QCD. Three of them have masses below the intrinsic 
scale of QCD Aqcd, and are called light. The other three have masses larger than Aqcd and are 
called heavy. Therefore, rather than describing and classifying EFT in general we will describe heavy 
quark systems and the EFT that can be constructed for them (as it will be more adequate for our 
purposes here). 

1.3 Heavy quark and quarkonium systems 

Three of the six quarks present in QCD have masses larger than Aqcd and are called heavy quarks. 
The three heavy quarks are the charm quark, the bottom quark and the top quark. The EFT will 
take advantage of this large mass of the quarks and construct an expansion in the heavy quark 
limit, of infinite quark masses. The simpler systems that can be constructed involving heavy quarks 
are hadrons composed of one heavy quark and one light (anti-)quark. The suitable effective theory 
for describing this kind of systems is known as Heavy Quark Effective Theory (HQET), and it is 
nowadays (together with chiral perturbation theory, which describes low energy interactions among 
pions and kaons, and the Fermi theory of weak interactions, which describes weak disintegrations 
below the mass of the W) a widely used example to show how EFT work in a realistic case (see £Q 
for a review of HQET). In brief, the relevant scales for this kind of systems are the heavy quark 
mass m and Aqcd- The effective theory is then constructed as an expansion in Aqcd/tti- The 
momentum of a heavy quark is decomposed as 

p = mv + k (1-2) 

where v is the velocity of the hadron (which is basically the velocity of the heavy quark) and k is 
a residual momentum of order Aqc d ■ The dependence on the large scale m is extracted from the 
fields, according to 

Q(x) = e- im ^ x Q v {x) = e - im ^- x [h v (x) + H v (x)} (1.3) 

and a theory for the soft fluctuations around the heavy quark mass is constructed. The leading 
order Lagrangian of HQET is given by 

C-hqet = h v iv ■ Dh v (1.4) 

This leading order Lagrangian presents flavor and spin symmetries, which can be exploited for 
phenomenology. 

The systems in which this thesis mainly focus (although not exclusively) are those known as 
heavy quarkonium. Heavy quarkonium is a bound state composed of a heavy quark and a heavy 
antiquark. We can therefore have charmonium (cc) and bottomonium (bb) systems. The heaviest 
of the quarks, the top, decays weakly before forming a bound state; nevertheless t — t production 
in the non-relativistic regime (that is near threshold) can also be studied with the same techniques. 
The relevant physical scales for the heavy quarkonium systems are the heavy quark mass m, the 
typical three momentum of the bound state mv (where v is the typical relative velocity of the 
heavy quark-antiquark pair in the bound state) and the typical kinetic energy mv 2 . Apart from 
the intrinsic hadronic scale of QCD Aqcd- The presence of all those scales shows us that heavy 
quarkonium systems probe all the energy regimes of QCD. From the hard perturbative region to 
the low energy non-perturbative one. Heavy quarkonium systems are therefore an excellent place 
to improve our understanding of QCD and to study the interplay of the perturbative and the non- 
perturbative effects in QCD 2 . To achieve this goal, EFT for this system will be constructed. 
Using the fact that the mass m of the heavy quark is much larger than any other scale present 



18 



Chapter 1. General introduction 



in the problem (a procedure which is referred to as integrating out the scale m) one arrives at an 
effective theory known as Non-Relativistic QCD (NRQCD) 3 . In that theory, which describes the 
dynamics of heavy quark-antiquark pairs at energy scales much smaller than their masses, the heavy 
quark (and antiquark) is treated non-relativistically by (2 component) Pauli spinors. Also gluons 
and light quarks with a four momentum of order m are integrated out and not present any more 
in the effective theory. What we have achieved with the construction of this EFT is the systematic 
factorization of the effects at the hard scale m from the effects coming from the rest of scales. 
NRQCD provides us with a rigorous framework to study spectroscopy, decay, production and many 
other heavy quarkonium processes. The leading order Lagrangian for this theory is given by 



where tp is the field that annihilates a heavy quark and x the field that creates a heavy antiquark. 
Sub- leading terms (in the 1/m expansion) can then be derived. One might be surprised, at first, 
that heavy quarkonium decay processes can be studied in NRQCD. Since the annihilation of a QQ 
pair will produce gluons and light quarks with energies of order m, and those degrees of freedom are 
not present in NRQCD. Nevertheless those annihilation processes can be explained within NRQCD 
(in fact the theory is constructed to reproduce that kind of physics) . The answer is that annihilation 
processes are incorporated in NRQCD through local four fermion operators. The QQ annihilation 
rate is represented in NRQCD by the imaginary parts of QQ — > QQ scattering amplitudes. The 
coefficients of the four fermion operators in the NRQCD Lagrangian, therefore, have imaginary 
parts, which reproduces the QQ annihilation rates. In that way we can describe inclusive heavy 
quarkonium decay widths to light particles. 

NRQCD has factorized the effects at the hard scale m from the rest of scales in the problem. 
But if we want to describe heavy quarkonium physics at the scale of the binding energy, we will face 
with the complication that the soft, mv, and ultrasoft, mv 2 , scales are still entangled in NRQCD. 
It would be desirable to disentangle the effects of those two scales. To solve this problem one can 
proceed in more than one way. One possibility is to introduce separate fields for the soft and ultrasoft 
degrees of freedom at the NRQCD level. This would lead us to the formalism now known as velocity 
NRQCD (vNRQCD) g). Another possibility is to exploit further the non-relativistic hierarchy of 
scales in the system (m 3> mv 3> mv 2 ) and construct a new effective theory which just contains 
the relevant degrees of freedom to describe heavy quarkonium physics at the scale of the binding 
energy. That procedure lead us to the formalism known as potential NRQCD (pNRQCD) 5 and is 
the approach that we will take here, in this thesis. When going from NRQCD to pNRQCD one is 
integrating out gluons and light quarks with energies of order of the soft scale mv and heavy quarks 
with energy fluctuations at this soft scale. This procedure is sometimes referred to as integrating 
out the soft scale, although the scale mv is still active in the three momentum of the heavy quarks. 
The resulting effective theory, pNRQCD, is non-local in space (since the gluons are massless and the 
typical momentum transfer is at the soft scale) . The usual potentials in quantum mechanics appear 
as Wilson coefficients of the effective theory. This effective theory will be described in some more 
detail in section |3~T1 

The correct treatment of some heavy quark and quarkonium processes will require additional 
degrees of freedom, apart from those of HQET or NRQCD. When wc want to describe regions of 
phase space where the decay products have large energy, or exclusive decays of heavy particles, for 
example, collinear degrees of freedom would need to be present in the theory. The interaction of 
collinear and soft degrees of freedom has been implemented in an EFT framework in what now is 
known as Soft- Collinear Effective Theory (SCET) j5J |7] . This effective theory will also be described 
in a following section I5~2l Just let us mention here that, due to the peculiar nature of light cone 
interactions, this EFT will be non-local in a light cone direction (collinear gluons can not interact 




(1.5) 
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with soft fermions without taking them far off-shell). 

The study of heavy quark and quarkonium systems has thus lead us to the construction of 
effective quantum field theories of increasing richness and complexity. The full power of the quantum 
field theory techniques (loop effects, matching procedures, resummation of large logarithms...) is 
exploited to obtain systematic improvements in our understanding of those systems. 

1.4 Structure of the thesis 

This thesis is structured in the following manner. Next chapter (chapter |2J) is a summary of the 
whole thesis written in Catalan (it does not contain any information which is not present in other 
chapters, except for the translation). Chapter [3] contains an introduction to potential NRQCD and 
Soft-Collinear Effective Theory, the two effective theories that are mainly used during the thesis. 
The three following chapters (chapters El GO and HJ) comprise the original contributions of this thesis. 
Chapter0]is devoted to the study of the (infrared dependence of the) QCD static potential, employing 
pNRQCD techniques. ChapterElis devoted to the calculation of an anomalous dimension in SCET 
(two loop nf terms are obtained), which is relevant in many processes under recent study. And 
chapter [fj] is devoted to the study of the semi-inclusive radiative decays of heavy quarkonium to 
light hadrons, employing a combination of pNRQCD and SCET. Chapter is devoted to the final 
conclusions. This chapter is followed by three appendices. The first appendix contains definitions 
of several factors appearing throughout the thesis. The second appendix contains Feynman rules 
for pNRQCD an SCET. And, finally, the third appendix contains the factorization formulas for the 
NRQCD matrix elements in the strong coupling regime. 



Chapter 2 

Summary in Catalan 



Per facilitar la lectura, i una eventual comparacio amb d'altres referencies escrites en angles, incloem, 
en la taula l2~Tl la traduccio emprada per alguns dels termes presents en la tesi. 

2.1 Introduccio general 

Aquesta tesi versa sobre l'estudi de l'estructura i les interaccions dels constituents fonamentals de la 
materia. Varem arribar al final del segle XX descrivint les propietats mes fonamentals conegudes de 
la materia en termes de teories quantiques de camps (pel que fa a les interaccions electromagnetiques, 
nuclear forta i nuclear feble) i de la relativitat general (pel que fa a la interaccio gravitatoria) . El 
Model Estandard (ME) de les interaccions fonamentals en la natura engloba teories quantiques 
de camps per descriure les interaccions electromagnetiques (l'anomenda ElectroDinamica Quantica, 
EDQ) i nuclears febles (que estan unificades en l'anomenada teoria electro-feble) i per descriure les 
interaccions fortes (l'anomedada CromoDinamica Quantica, CDQ). Aquest ME ve complementat per 
la teoria classica (no quantica) de la gravitacio, la relativitat general. Tots els experiments que s'han 
dut a terme en acceleradors de partfcules (per tal d'estudiar els constituents basics de la materia), 
fins el dia d'avui, son consistents amb aquest marc teoric. Aixo ens ha portat a comengar el segle 
XXI esperant que la segiient generacio d'experiments destapi la ffsica que hi pot haver mes enlla 
d'aquestes teories, que han quedat ja ben establertes. Hi ha grans esperances posades en el gran 
accelerador hadronic, anomenat Large Hadron Collider (LHC), que s'esta construint actualment al 
CERN. Esta planificat que aquesta maquina comenci a ser operativa l'any 2007. El que s'espera es 
que l'LHC ens obri el camf cap a nous fenomens fisics no observats fins ara. Aixo es pot aconseguir 
de dues maneres. Una possibilitat es que l'LHC descobreixi la parti'cula de Higgs (l'unica parti'cula 
del ME que encara no s'ha observat) i que aixo desencadeni la descoberta de noves partfcules mes 
enlla del ME. L'altra possibilitat es que l'LHC demostri que no hi ha tal parti'cula de Higgs; cosa que 
demanaria un marc teoric totalment nou i diferent l'actual (per explicar les interaccions fonamentals 
de la natura) 1 . Les possibles extensions del ME han estat ja estudiades ampliament i amb gran detail. 
El que s'espera es que tots aquests efectes es facin palesos en aquesta nova generacio d'experiments. 
No cal dir que els experiments basats en acceleradors de partfcules no son Tunica opcio que tenim, 
per tal de descobrir efectes associats a nova ffsica. Una altra gran oportunitat (que tambe ha estat 
ampliament estudiada) es la d'observar el eel i la informacio que ens arriba d'ell (partfcules altament 
energetiques, fons cosmics de radiacio...). 

1 Intentarem no preocupar-nos gaire per la possibilitat que l'LHC descobreixi el Higgs, tanqui el ME i mostri que 
no hi ha efectes de nova ffsica en cap escala d'energia que serem capacos d'assolir amb acceleradors de partfcules. Tot 
i que aixo es possible no es, obviament, gens desitjable. 
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Pero en aquest canri a la recerca de la segiient teoria mes fonamcntal concguda fins al moment, no 
volcm perdre Fhabilitat de fer sevir aquesta teoria per fer prediccions concises per un ampli ventall 
de processos fi'sics, i tambe volem poder entendre (d'una forma no ambigua) com la teoria precedent 
es pot obtenir a partir de la nova. Obviament el somni de qualsevol fi'sic es trobar una descripicio 
unificada de les quatre interaccions fonamentals conegudes de la materia; pero no al preu de tenir 
una teoria que pot explicar-ho tot pero que es tant complicada que, de fet, no explica res. Acabarem 
aquests paragrafs que fan de prefaci a la tesi amb un petit acudit. D'acord amb el que hem dit aquf, 
el titol de la conferencia que suposaria el punt i final de la fisica (teorica) no es: La teoria M: una 
unificacio de totes les interaccions de la natura, sino mes aviat: Com obtenir el metabolisme d'una 
vaca a partir de la teoria M. 

2.1.1 Teories Efectives 

En aquesta tesi ens centrarem en l'estudi de sistemes que involucren el sector de les interaccions 
fortes en el ME. La part del ME que descriu les interaccions fortes es, com s'ha comentat abans, 
la CromoDinamica Quantica. CDQ es una teorica quantica de camps basada en el grup no abelia 
SU (3) i descriu les interaccions entre quarks i gluons. El seu Lagrangia es extremadament simple i 
ve donat per 

Nf 

Cqcd = ]T ft (ip- rat) q t -G^ a G; u (2.1) 

i=l 

En aquesta equacio qi son els camps associats als quarks, igG^ v = [D^, D v \, amb = d^+igA^, 
son els camps pels gluons i Nf es el niimero total de sabors (tipus) de quarks. La CDQ presenta les 
propietats de Uibertat asimptotica i de confinament. La constant d'acoblament de les interaccions 
fortes esdeve gran a energies petites i tendeix a zero per energies grans. D'aquesta manera, per 
energies altes els quarks i els gluons es comporten com a particules lliurcs, mentre que a baixes 
energies apareixen sempre confinats a l'interior d'hadrons (en una combinacio singlet de color). La 
CDQ desenvolupa una escala intrmseca, Aqcd, a baixes energies; escala que dona la contribucio 
principal a la massa de la majoria dels hadrons. Aqcd es pot interpretar de diferents maneres, pero 
es basicament l'escala d'energia on la constant d'acoblament de les interaccions fortes esdeve d'ordre 
1 (i la teoria de perturbacions en a s ja no es fiable). Es pot pensar que es una escala de l'ordre 
de la massa del proto. La presencia d'aquesta escala intrinseca i el fet, intimament relacionat, que 
l'espectre de la teoria consisteixi en estats hadronics singlets de color (i no dels quarks i gluons) 
provoca que els calculs directes des de CDQ siguin extremadament complicats, sino impossibles, per 
molts sistemes fisics d'interes. Les tecniques conegudes amb el nom de teories efectives (TE) ens 
ajudaran en aquesta tasca. 

Com a regla general, l'estudi de qualsevol proces, en teoria quantica de camps, que involucri 
mes d'una escala fisica rellevant es complicat. Els calculs (i les integrals) que ens apareixeran poden 
resultar molt complicats si mes d'una escala entra en ells. La idea sera doncs construir una nova teoria 
(la teoria efectiva), derivada de la teoria fonamental, de manera que nomes involucri els graus de 
Uibertat rellevants per la regio que ens interessa. La idea general que hi ha sota les tecniques de TE es 
simplement la segiient: per tal d'estudiar la fisica d'una determinada regio d'energies no necessitem 
coneixer la dinamica de les altres regions de forma dctallada. Aquest es, obviament, un fet ben 
conegut i ampliament acceptat. Per exemple, tothom entent que per descriure una reaccio qufmica 
no cal coneixer la interaccio quantica electrodinamica entre els fotons i els electrons, per contra un 
model de l'atom que consisteixi en un nucli i electrons orbitant al voltant es mes convenient. I de la 
mateixa manera no cal usar aquest model atomic per tal de descriure un proces biologic macroscopic. 
La implementacio d'aquesta ben coneguda idea en el marc de la teoria quantica de camps es el que 
es coneix sota el nom generic de Teories Efectives. Tal i com ja s'ha dit abans, aquestes tecniques 
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esdevenen especialment utils en l'estudi de processes que involucren les interaccions fortes. Per tal de 
construir una TE cal seguir els segiients passos generals (a grans trets) . En primer Hoc es necessari 
identificar els graus de llibertat que son rellevants pel problema en que estem interessats. Despres 
cal for us de les simetries presents en el problema i, finalment, hem d'aprofitar qualsevol jerarquia 
d'escales que hi pugui haver. Es important remarcar que el que estem font no es construir un model 
pel proces que volem estudiar. Per contra la TE esta construida de manera de sigui equivalent a 
la teoria fonamental, en la regio on es valida; estem obtenint els resultats desitjats a partir d'una 
expansio ben controlada de la nostra teoria fonamental. 

Mes concretament, en aquesta tesi ens centrarem en l'estudi de sistemes que involucren els 
anomenats quarks pesats. Com es ben conegut hi ha sis sabors (tipus) de quarks en CDQ. Tres 
d'ells tene masses per sota de l'escala Aqcd i s'anomenen lleugers, mentre que els altres tres tenen 
masses per sobre d'aquesta escala Aqcd i s'anomenen pesats. El que farem a continuacio es descriure 
sistemes amb quarks pesats i les teories efectives que es poden construir per ells. 

2.1.2 Sistemes de quarks pesats i quarkoni 

El que faran les TE pels sistemes amb quarks pesats es aprofitar-se d'aquesta escala gran, la massa, 
i construir una expansio en el limit de quarks infim'tament massius. Els sistemes mes simples que 
es poden tenir involucrant quarks pesats son aquells composats d'un quark pesat i un (anti-)quark 
llcuger. La TE adequada per descriure aquest tipus de sistemes rep el nom de Teoria Efectiva per 
Quarks Pesats (TEQP). Aquesta teoria es avui en dia, i juntament amb la teoria de perturbacions 
quiral (que descriu les interaccions de baixa energia entre pions i kaons) i la teoria de Fermi per 
les interaccions febles (que decriu les desintegracions febles per a energies per sota de la massa del 
boso W), un exemple ampliament usat per mostrar com les TE funcionen en un cas realista. De 
manera molt breu, les escales fisiques rellevants per aquest sistema son la massa m del quark pesat 
i Aqcd- La TE es construeix, per tant, com una expansio en Aqcd/™- El moment del quak pesat 
es descomposa d'acord amb 

p = mv + k (2.2) 

on v es la velocitat de l'hadro (que es basicament la velocitat del quark pesat) i k es un moment 
residual d'ordre Aqcd- La dependencia en l'escala m s'extreu dels camps de la TE d'acord amb 

Q{x) = e- im ^ v - x Q v {x) = e - im ^ x [h v (x) + H v (x)} (2.3) 

i es construeix una teoria per les fluctuations suaus al voltant de la massa del quark pesat. El 
Lagrangia de la TEQP a ordre dominant ve donat per 

Chqet = Kiv ■ Dh v (2.4) 

Aquest Lagrangia presenta simetries de sabor i spin, que es poden aprofitar per a fer fenomenologia. 

Els sistemes en que aquesta tesi se centrara (encara que no de manera exclusiva) son aquells 
coneguts amb el nom de quarkoni pesat. El quarkoni pesat es un estat lligat composat per un quark 
pesat i un antiquark pesat. Per tant podem tenir sistemes de charmoni (cc) i de bottomoni (bb). El 
mes pesat de tots els quarks, el quark top, es dcsintcgra a traves de les interaccions febles abans que 
pugui formar estats Uigats; de tota manera la produccio de parelles t — t prop del llindar de produccio 
(per tant, en un regim no relativista) es pot estudiar amb les mateixes tecniques. Les escales fisiques 
rellevants pels sistemes de quarkoni pesat son l'escala m de la massa del quark pesat, el tri-moment 
tipic de l'estat lligat mv (v es la velocitat relativa tipica de la parella quark-antiquark en l'estat 
lligat) i l'energia cinetica tipica mv 2 . A part de l'escala intrinseca de la CDQ, Aqcd, que sempre 
es present. La presencia simultania de totes aquestes escales ens indica que els sistemes de quarkoni 
pesat involucren tots els rangs d'energia de CDQ, des de les regions perturbatives d'alta energia fins 
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a les no-perturbatives de baixa energia. Es per tant un bon sistema per estudiar la interaccio entre 
els efectes perturbatius i els no perturbatius en CDQ i per millorar el nostre coneixement de CDQ 
en general. Per tal d'aconseguir aquest objectiu contruirem TE adequades per la descripcio d'aquest 
sistema. Si fern servir el fet que la massa m es molt mes gran que qualsevol altra escala d'energia 
present el problema, arribem a una TE coneguda amb el nom de CDQ No Relativista (CDQNR). En 
aquesta teoria, que descriu la dinamica de parelles de quark-antiquark per energies forga menors a les 
seves masses, els quarks pesats venen representats per spinors no relativistes de dues components. A 
mes a mes, gluons i quarks lleugers amb quadri-moment a l'escala m son intergats de la teoria i ja no 
hi apareixen. El que hem aconseguit amb la construccio d'aquesta teoria es factoritzar, de manera 
sistematica, els efectes que venen de l'escala m de la resta d'efectes provinents de les altres escales del 
problema. CDQNR ens proporciona un marc teoric rigoros on estudiar processos de desintegracio, 
produccio i espectroscopia de quarkoni pesat. El Lagrangia a ordre dominant ve donat per 

Cnrqcd = V f (iD + ^D 2 ) V + X f (iD - ^D 2 ) X (2.5) 

on %j) es el camp que anihila el quark pesat i \ e l camp que crea l'antiquark pesat. Termes sub- 
dominants, en l'expansio en 1/m, poden ser derivats. D'entrada pot resultar sorprenent que els 
processos de desintegracio puguin ser estudiats en el marc de la CDQNR. L'anihilacio de la parella 
QQ produira gluons i quarks lleugers amb energies d'ordre m, i aquests graus de llibertat ja no son 
presents en CDQNR. Tot i aixo els processos de desintergacio poden ser estudiats en el marc de la 
CDQNR, de fet la teoria esta construida per tal de podcr explicar aquests processos. La resposta 
es que els processos d'anihilacio s'incorporen en CDQNR a traves d'interaccions locals de quatre 
fermions. Les raons de desintegracio venen representades en CDQNR per les parts imaginaries de 
les amplituds de dispersio QQ — ► QQ. Els coeficients dels operadors de quatre fermions tenen, per 
tant, parts imaginaries que codifiquen les raons de desintegracio. D'aquesta manera podem estudiar 
les desintegracions inclusives de quakonium pesat en partfcules lleugeres. 

La CDQNR ens ha factoritzat els efectes a l'escala m de la resta. Ara be, si volem estudiar la fisica 
del quarkoni pesat a l'escala de l'energia de lligam del sistema, ens trobarem amb el problema que 
les escales suau, corresponent al tri- moment tipic mv, i ultrasuau, corresponent a l'energia cinetica 
tipica mv 2 , estan encara entrellagades en CDQNR. Seria desitjable separar els efectes d'aquestes dues 
escales. Per tal de solucionar aquest problema es pot procedir de mes d'una manera. L'estrategia 
que emprarem en aquesta tesi es la d'aprofitar de manera mes amplia la jerarquia no relativista 
d'escales que presenta el sistema (m 3> mv 3> mv 2 ) i construir una nova teoria efectiva que nomcs 
contingui els graus de llibertat rellevants per tal de descriure els sistemes de quarkoni pesat a 
l'escala de l'energia de lligam. La teoria que s'obte es coneguda amb el nom de CDQNR de potencial 
(CDQNRp). Aquesta teoria sera descrita breument en la segiient seccio. 

El tractament correcte d'alguns sistemes de quarks pesats i de quarkoni pesat demanara la 
presencia de graus de llibertat addicionals, a part dels presents en TEQP o en CDQNR. Quan volem 
descriure regions de l'espai fasic on els productes de la desintegracio tenen una energia gran, o quan 
volguem descriure desintegracions exclusives, per exemple, graus de llibertat col-lineals hauran de 
ser presents en la teoria. La interaccio dels graus de llibertat col-lineals amb els graus de llibertat 
suaus ha estat implementada en el marc de les TE en el que avui es coneix com a Teoria Efectiva 
Col-lineal-Suau (TECS). Aquesta teoria la descriurem tambe breument en la segiient seccio. 

En definitiva, l'estudi de sistemes de quarks pesats i quarkoni ens ha portat a la construccio 
de teories efectives de camps de riquesa i complexitat creixents. Tota la potcncia de les tccniqucs 
de teoria quantica de camps (efectes de bagues, resumacio de logartimes...) es explotat per tal de 
millorar la nostra comprensio d'aquests sistemes. 
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Angles 


Oatala 


Quantum Chrome-Dynamics (QCD) 


CromoDinamica Quantica (CDQ) 


bolt-Gollmcar riliective Ineory (aGEl) 


leoria lelectiva Col-lmeal-buau (llLGb) 


loop 


baga 


btandard Model (bM) 


Model iestandard (IVLfcj) 


Quantum ElectroDynamics (QED) 


Electrodinamica Quantica (EDQ) 


quarkonmm 


quarkom 


Heavy (^uark reflective Ineory (H(4-Cj1 j 


leoria rjicctiva per (Quarks Fesats (l-b^r) 


Non-Rclativistic QCD (NRQCD) 


CDQ No Rclativista (CDQNR) 


potential NRQCD (pNRQCD) 


CDQNR de potencial (CDQNRp) 


matching coefficients 


coeficients de coincidencia 


label operators 


operadors etiqueta 


jet 


doll 



Taula 2.1: Traduccio anglcs-catala d'alguns termes usats en la tesi. 

2.2 Rerefons 
2.2.1 CDQNRp 

Com ja s'ha dit abans, les escales rellevants pels sistemes de quarkoni pesat son la massa m, l'escala 
suau mv i l'escala ultrasuau mv 2 . A part de l'escala Aqcd- Quan aprofitcm la jcrarquia no rclativista 
del sistema en la seva totalitat arribem a la CDQNRp. Per tal d'identificar els graus de llibertat 
rellevants en la teoria final, cal especificar la importancia relativa de Aqcd respecte les escales suau i 
ultrasuau. Dos regims rellevants han estat identificats. Son els anomenats regim d'acoblament feble 
mv 2 > Aqcd i regim d'acoblament fort mv > Aqcd 3> mv 2 . 

Regim d'acoblament feble 

En aqucst regim els graus de llibertat de CDQNRp son semblants als de CDQNR, pero amb les 
cotes superiors en energia i tri-moments abaixades. Els graus de llibertat de CDQNRp consisteixen 
en quarks i antiquarks pesats amb un tri-moment fitat superiorment per v p (|p| <C v p <C m) i una 

energia fitada per v us <C v us <C |p|), i en gluons i quarks lleugers amb un quadri- moment fitat 
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per v us . El Lagrangia es pot escriure com 

CpNRQCD = J d 3 r Tr |s f (id Q - h,(r, p, P R , S u S 2 )) S + f (i£) - fc (r, p, P R , Si, S 2 )) o| + 

+V A (r)Tr {0 f r • 5 E S + S f r • 5 E 0} + ^^Tr {O f r • . 9 E O + f 0r • gE} - 



amb 



-I G ^ a + f>^ (2.6) 



i=l 



/* s (r, p, P R , Si, S 2 ) = + + ^( r , p, P R , Si, S 2 ) (2.7) 

m rod 2 mtot 

/»„(r, p, P R , Si, S 2 ) = -H- + - + V (r, p, P R , Si, S 2 ) (2.8) 
m ro d 2 m to t 

i 

L» O = id O - g[A (R, t), O] P R = -iD R m rcd = — — - m to t = rn x + m 2 (2.9) 

m to t 

S es el camp singlet pel quarkoni i O el camp octet per ell. E representa el camp cromoelectric. Po- 
dem veure que els potencials usuals de mecanica quantica apareixen com a coeficients de coincidencia 
en la teoria efcctiva. 



Regim d'acoblament fort 

En aquesta situacio la fisica a l'escala de l'energia de lligam esta per sota de l'escala Aqcd- Per 
tant es millor discutir la teoria en tcrmes de graus de llibertat hadronics. Guiant-nos per alguncs 
consideracions generals i per indicacions provincnts CDQ en el reticle, podem suposar que el quarkoni 
ve descrit per un camp singlet. I si ignorem els bosons de Goldstone, aquests son tots els graus de 
llibertat en aques regim. El Lagrangia ve ara donat per 

^pNrqcd — J d 3 R J d 3 r St(ia D _fe a ( Xl)X2> p lj p 2j s 1 ,S 2 ))S (2.10) 

amb 

P 2 P 2 

/i s (xi,x 2 ,pi,p 2 ,Si,S 2 ) = — — + —2- + V s (xi,x 2 ,pi,p 2 ,Si,S 2 ) (2.11) 

El potencial V s es ara una quantitat no perturbativa. El procediment de coincidencia de la teoria 
fonamental i la teoria efectiva ens donara expressions pel potencial (en termes de les anomenades 
bagues de Wilson). 



2.2.2 TECS 

L'objectiu d'aquesta teoria es descriure processos on graus de llibertat molt energetics (col-lineals) 
interactuen amb graus de llibertat suaus. Aixi la teoria es pot aplicar a un ampli vcntall de processos, 
on aquesta situacio cinematica es present. Qualsevol proces que contingui hadrons molt energetics, 
juntament amb una font per ells, contindra particules, anomenades col-lineals, que es mouen a 
prop d'una direccio del con de Hum n^. Com que aquestes particules han de tenir una energia E 
gran i alhora una massa invariant petita, el tamany de les components del seu quadri-moment (en 
coordenades del con de Hum, = (hp)n^/2 +p' J L + (np)n M /2) es molt diferent. Tipicament np <~ E, 
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p± ~ EX i np <~ _EA 2 , amb A un parametre petit. Es d'aquesta jerarquia d'escales que la TE treura 
profit. 

Els graus de llibertat que cal inlcoure en la teoria efectiva depenen de si un vol estudiar processos 
inclusius o exclusius. Les dues teories que en resulten es coneixen amb els noms de TECSi i TECSn, 
respectivament. 



TECSi 

Aquesta es la teoria que conte graus de llibertat col-lineals (p M = (np,p± 7 np) ~ (1,A, A 2 )) i ultra- 
suaus (p p ~ (A 2 , A 2 , A 2 )). Els graus de llibertat col-lineals tenen massa invariant d'ordre EAqcd- 
Malauradament, tant en TECSi com en TECSn, no hi ha una notacio estandard en la litcratura. 
Dos formalismes (suposadament equivalents) han estat usats. 
El Lagrangia a ordre dominant ve donat per 

£c = W {inD + gnA n , q + (f ± + g^ q ) wlw^ (f ± + gjfc^) } |^ n , p (2.12) 

£n tP es el camp pel quark col-lineal, A nyP el camp pel gluo col-lineal (la dependencia en les escales 
grans ha estat extreta d'ells de manera scmblant a en TEQP). Els V son els anomenats operdors 
ctiqueta que donen les components grans (extretes) dels camps. Les W son lmies de Wilson. 



TECS„ 

Aquesta es la teoria que descriu processos on els graus de llibertat col-lineals en l'estat final tenen 
massa invariant d'ordre Aq C£) . Aquesta teoria es mes complicada que l'anterior, ja que en el proces 
d'anar des de CDQ a TECSn la presencia de dos tipus de modes col-lineals s'ha de tenir en comptc. 
En la tesi basicament no usarem aquesta teoria i, per tant, no en direm res mes. 



2.3 El potencial estatic singlet de CDQ 

El potencial estatic entre un quark i un antiquark es un objecte clau per tal d'entendre la dinamica 
de la CDQ. Aqui ens centrarem en estudiar la dependencia infraroja del potencial estatic singlet. 
Obtindrcm la dependencia infraroja sub-dominant del mateix fent servir la CDQNRp 
L'expansio perturbativa del potencial estatic singlet ve donada per 



Vj°\r) = - 



C f a s (l/r) 



a g (l/r) 

47T 



(ai + 2 7B /3o) + 



47T 



a 2 - 



+ -5" + 4 Te )ft > + IE (4oiA, + 2/3i) + 



on 



47T 



«3 



a 1 = —C A - Y Tpn f 



C\\ogrn 



(2.13) 
(2.14) 



a-2 = 



4343 
162 



4tt 2 - 



22 



C(3) 



Ci 



1798 56 



81 



f - 16C(3) 



CfT F n f + [ —T F n f 



20, 



C(3) 



CaTfuj 



(2.15) 
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el logaritme que veiem en l'expressio pel potencial es la dependencia infraroja dominant. Aqui 
trobarem la dependencia infraroja sub-dominant; es a dir una part de la correccio a quart ordre del 
potencial. Per fer-ho estudiarem el proces de fer coincidir CDQNR amb CDQNRp. El que cal fer es 
calcular la conicidencia a ordre r 2 en l'expansio multipolar de CDQNRp. Per fer aixo cal evaluar el 
segon diagrama de la part dreta de la igualtat de la figura^21 Quan calculem la primera correccio en 
a s d'aquest diagrama (despres del terme dominant) obtenim la dependencia infraroja sub-dominant 
que busquem (el terme dominant del diagrama donava la dependencia infraroja dominant). El 
resultat pels termes infrarojos sub-dominants del potencial es 

Vj 0) (r) = CEa TFm - 
C f a s (l/r) ^ s (l/r) V 16tt 2 / 11 2 \ 2 

1 1 ™ A [-—C A + -n f \\ogrii- (2.16) 



4tt/ 3 ^ V 3 3 



C f a s {l/r) f a s (l/r) Y I6n 

47T 



—^-C A (at + 2 1E /3 Q - — (20ti/ - Ca(12tt 2 + 149)) J logr/a (2.17) 



2.4 Dimensio anomala del corrent lleuger-a-pesat en TECS 
a dues bagues: termes nj 

Els corrents hadronics lleuger-a-pesat Jhad = qTb (b representa el quark pesat i q el quark lleuger), 
que apareixen en operadors de la teoria nuclear feble a una escala d'energia p. ~ m&, es poden fer 
coincidir amb els corrents de TECSi. A ordre mes baix en el parametre d'expansio A el corrent en 
TECS ve donat per 

J hSi ET = c o (np, n) £n, P Th + a (np, nqi,fx) £„, p {gnA nm )Th H (2-18) 

Es a dir, un nombre arbitrari de gluons nA n ^ q poden ser afegits, sense que aixo suposi supressio en el 
comptatge en el parametre A. Els coeficients de Wilson poden ser evolucionats, en la teoria efectiva, 
a una escala d'energia mes baixa. Com que tots els corrents estan relacionats per invariancia de 
galga col-lineal, es suficient estudiar el corrent ^T/i (que es obviament mes simple). L'evolucio del 
corrent a una baga ve determinada per la dimensio anomala 




P es el moment total sortint del doll de partfcules. Aqui volem trobar els termes n/ de la correccio 
a dues bagues d'aquest resultat. Per tal de calcular-los cal evaluar els diagrames de la figura l5~31 A 
part tambe necessitem la correccio a dues bagues dels propagadors del quark col-lineal i del quark 
pesat. La correccio del propagador del quark col- lineal coincideix amb la usual de CDQ (ja que 
en el seu calcul nomes hi entren partfcules col-lineals, i no ultrasuaus); mentre que la correccio al 
propagador del quark pesat es la usual de TEQP. Tenint en compte el resultat dels diagrames i 
aquestes correccions als propagadors, obtenim el resultat desitjat pels termes n/ a dues bagues de 
la dimensio anomala 



<2 ba9 ues n f) = ( ^ ^ F U fCf ( 125 ^ 20 / \ 

' Utt/ 3 V 18 2 3 ^nPJ J 



(2.20) 
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2.5 Desintegracions radiatives de quarkoni pesat 



Les desintegracions semi-inclusives radiatives de quarkoni pesat a hadrons lleugers han estat estu- 
diades des dels inicis de la CDQ. Aquests primers treballs tractaven el quarkoni pesat en analogia 
amb la desintegracio de l'orto-positroni en EDQ. Diversos experiments, que es van fer posteriormcnt, 
van mostrar que la regio superior z — > 1 de l'espectre del foto {z es la fraccio d'energia del foto, 
respecte la maxima possible) no podia ser ben explicada amb aquests calculs. Posteriors calculs 
de correccions relativistes i de resumacio de logaritmes, tot i que anaven en la bona direccio, no 
eren tampoc suficients per explicar les dades cxpcrimcntals. Per contra, l'espectre podia ser ben 
explicat amb models que incorporaven una massa pel gluo. L'aparicio de la CDQNR va permetre 
analitzar aquestes desintegracions de manera sistematica, pero, tot i aixf, una massa finita pel gluo 
semblava necessaria. Ben aviat, per aixo, es va notar que en aquesta regio superior la factoritzacio 
de la CDQNR no funcionava. S'havien d'introduir les anomenades funcions d'estrucutra (en el canal 
octet de color), que integraven contribucions de diversos ordres en Pexpansio de CDQNR. Alguns 
primers intents de modelitzar aquestes funcions d'estructura dugueren a resultats en fort desacord 
amb les dades. Mes endavant es va reconeixer que per tractar correctament aquesta regio superior 
de l'espectre calia combinar la CDQNR amb la TECS (ja que els graus de llibertat col-lineals tambe 
eren importants en aquesta regio cinematica). D'aquesta manera les resumacions de logaritmes van 
ser estudiades en aquest marc (i es corregiren i ampliaren els resultats previs). Aqui farem servir 
una combinacio de la CDQNRp amb la TECS per tal de calcular aquestes funcions d'estructura 
suposant que el quarkoni que es desintegra es pot tractar en el regim d'acoblamcnt feble. Quan 
combinem de manera consistent aquests resultats amb els resultats previs coneguts, s'obte una bona 
descripcio de l'espectre (sense que ja no calgui introduir una massa pel gluo) en tot el rang de z. 

Per tal de calcular aquestes funcions d'estrucutra, el primer que cal fer es escriure els corrents 
en CDQNRp+TECS, que es on els calcularem. Un cop es te aixo ja es poden calcular els diagrames 
corresponents i aleshores, comparant amb les formules de factoritzacio per aquest proces, es poden 
indentificar les funcions d'estrucutra desitjades. Els diagrames que cal calcular venen representats 
a la figura RT2l Del calcul d'aquests diagrames s'obtenen les funcions d'estructura. El resultat que 
s'obte es divergent ultraviolat i ha de ser renormalitzat. Un cop s'ha fet aixo, si comparem el resultat 
teoric que tenim ara per l'espectre amb les dades experimentals en la regio superior, trobem un bon 
acord; tal i com es pot veure en la figura 16.51 (les dues corbes en la figura representen diferents 
esquemes de renormalitzacio). Fins ara hem pogut explicar, doncs, la regio superior de l'espectre. 
El que ara falta fer es veure si aquests resultats es poden combinar amb els calculs anteriors, per la 
resta de l'espectre, i obtenir un bon acord amb les dades experimentals en tot el rang de z. Cal anar 
amb compte a l'hora de combinar aquests resultats, ja que en les diferents regions de l'espectre son 
necessaries diferents aproximacions teoriques (per tal de poder calcular). El proces emprat consisteix 
doncs en expandir (per z en la regio central) les expressions que hem obtingut per la regio superior 
de l'espectre. Aleshores cal combinar les expressions d'acord amb la formula 
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(2.21) 



on SC representa la contribucio sinlget de color, OC la contribucio octet de color i els superfndexs c i e 
es refereixen a les expressions per la regio central i per l'extrem superior de l'espectre, respectivament. 
Quan usem aquesta formula aconseguim obtenir l'expressio valida per la regio central en la regio 
central i l'expressio valida per l'extrem superior de l'espectre en l'extrem superior, a part de termes 
que son d'ordre superior en el comptatge de la teoria efectiva en les respectives regions. I ho hem fet 
sense haver d'introduir talis o cotes arbitraries per tal de delimitar les diferents regions de l'espectre 
(cosa que hagues introduit incerteses teoriques basicament incontrolables en els nostres resultats). 
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Quan comparem el resultat d'aquesta corba 2 (que ara ja conte tots els termes que, d'acord amb el 
nostre comptatge, han de ser presents) amb les dades experimentals, obtenim un molt bon acord. 
La comparacio es pot veure a les figures 16.131 i lfi.141 (la corba vermella (clara) contmua en aquestes 
figures es la prediccio teorica per l'espectre) . 

Un cop ja tenim l'espectre ben descrit des del punt de vista teoric, podem fer-lo servir per 
cstudiar propietats del quarkoni pesat. En concret, es possible fer servir aquests espectres per tal 
de determinar en quin regim d'acoblament es troben els diferents quarkonis que es desintegren. Si 
calculem el quocient d'espectres de dos estats (n i r) en el regim d'acoblament fort obtenim 
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(totes les quantitats que apareixen en aquesta equacio son conegudes), mentre que si un dels dos 
estats es en el regim d'acoblament feble la formula que obtenim presenta una dependencia en z 
different. Per tant si la formula anterior reprodueix be el quocient d'espectres, aixo ens estara 
indicant que els dos quarkonis estan en el regim d'acoblament fort, mentre que si no es aixf almenys 
un dels dos sera en el regim d'acoblament feble. Com que hi ha dades disponiblcs pels estats 
T(15), T(2S) i T(3S) podem portar aquest proces a la practica. La comparacio amb els resultats 
experimentals es pot veure a les figures IB~T%1 16.191 i ET^UI (l'estat T(15) esperem que estigui en el 
regim d'acoblament feble, cosa que es compatible amb la grafica lfi.lSJl . Els errors son molt grans, 
pero T(2S) i T(3S) semblen compatibles amb ser estats de regim d'acoblament fort (cal comparar la 
corba contmua amb els punts. Si coincideixen indica que els dos estats son en el regim d'acoblament 
fort). 



2.6 Conclusions 

En aquesta tesi hem fet servir les tecniques de teories efectives per tal d'estudiar el sector de quarks 
pesats del Model Estandard. Ens hem centrat en l'estudi de tres temes. En primer Hoc hem estudiat 
el potencial estatic singlet de CDQ, fent servir la CDQ No Relativista de potencial. Amb l'ajuda 
d'aquesta teoria efectiva hem estat capacos de determinar la dependencia infraroja sub-dominant 
d'aquest potencial estatic. Entre altres possibles aplicacions, aquest resultat es rellevant en l'estudi 
de la produccio de t — t prop del llindar de produccio (a tercer ordre) . Aquest es un proces que 
cal ser estudiat amb molt de detail amb vista a la possible futura construccio d'un gran accelerador 
lineal electro-positro. Despres hem estudiat una dimensio anomala en la TECS. Aquesta teoria te 
aplicacions molt importants en el camp de la fi'sica de mesons B. I aquest es un camp de gran 
importancia per a la recerca indirecta de processos associats a nova fisica (mitjangant l'estudi de 
la violacio de CP i de la matriu de CKM). Finalment hem estudiat les desintegracions raditives 
semi-inlcusives de quarkoni pesat a hadrons lleugers. Per tal d'explicar be aquest proces ha estat 
necessaria una combinacio de la CDQNRp amb la TECS. Mirant-s'ho des de la perspectiva actual, 
es pot veure aquest proces com un bonic exemple de com, una vegada s'incorporen tots els graus de 
llibertat rellevants en un problema (i es fa servir un comptatge ben definit per ells), aquest es ben 
descrit per la teoria. Un cop aquest proces esta entes, es pot fer servir per estudiar algunes de les 
propietats del quarkoni pesat que es desintegra; com tambe hem mostrat en la tesi. 



2 Tambe cal afegir les anomcnades contributions de fragmentacio. A l'ordre en que estem treballant aquf son 
completament independents de les contributions directes de la formula anterior. 



Chapter 3 

Background 



In this chapter we describe the two effective field theories that will be mainly used and studied in 
the thesis: potential Non-Relativistic QCD and Soft-Collinear Effective Theory. It does not attempt 
to be a comprehensive review but just provide the sufficient ingredients to follow the subsequent 
chapters. 

3.1 potential Non Relativistic QCD 

As it has already been explained in the introduction of the thesis, heavy quarkonium systems are 
characterized by three intrinsic scales. Those are, the heavy quark mass m (which is referred to 
as the hard scale and sets the mass of the quarkonium state) , the relative three-momentum of the 
heavy quark-antiquark pair |p| ~ mv (which is referred to as the soft scale and sets the size of the 
bound state, v is the typical relative velocity between the quark and the antiquark) and the kinetic 
energy of the heavy quark and antiquark E ~ mv 2 (which is referred to as the ultrasoft scale and 
sets the binding energy of the quarkonium state), and by the generic hadronic scale of QCD Aqcd- 
All those scales are summarized in table 13.11 The interplay of Aqc d with the other three scales 
determines the nature of the different heavy quarkonium systems. By definition of heavy quark, 
m is always much larger than Aqcd] so the inequality m 3> Aqcd always holds. Exploiting the 
inequality m ^S> \p\,E one arrives at Non-Relativistic QCD (NRQCD), as it has been described in 
the previous chapter (note that at this level, after the definition of heavy quark, one still does not 
need to specify the interplay of Aqcd with the remaining scales, to identify the relevant degrees of 
freedom). Going one step further, using the full non-relativistic hierarchy of the heavy quarkonium 
systems m ^> mv ^> mv 2 , one arrives at potential NRQCD (pNQRCD) 1 . Now it is necessary to 
set the relative importance of Aqcd with the scales |p| and E to fix the degrees of freedom of the 
resulting theory, the aim of which is to study physics at the scale of the binding energy E. Two 
relevant regimes have been identified so far; the so called weak coupling regime, where mv 2 > Aqcd, 
and the so called strong coupling regime, where mv > Aqcd ^ mv 2 . 

3.1.1 Weak coupling regime 

In this situation, the degrees of freedom of pNRQCD are not very different from those of NRQCD. 
They are heavy quarks and antiquarks with a three momentum cut-off v p (|p| <C v p <C m) and an 

2 

energy cut-off v us (^ <C v us <C |p|), and gluons and light quarks with a four momentum cut-off v us . 

1 See [HI for a review of pNRQCD. 
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Scale 


Value 


hard 


m 


soft 


mv 


ultrasoft 


mv 2 


generic hadronic QCD scale 





Table 3.1: Relevant physical scales in heavy quarkonium 



The most distinct feature is that now non local terms in r, that is potentials, can appear (as it has 
been discussed before, in the introductory section fOjl . These degrees of freedom can be arranged 
in several ways in the effective theory. One first way is to express them with the same fields as in 
NRQCD. Then the pNRQCD Lagrangian has the following form 

LpNRQCD = VxRQCD + Lpot (3.1) 

where L^ R q CD is the NRQCD Lagrangian but restricted to ultrasoft gluons and L pot is given by 

Lpot = ~ J d 3 Xid 3 x 2 V' t {t, xi) x (t, x 2 ) V (r, pi, p 2 , Si, S 2 ) (its gluon fields)x f (t, x 2 ) ip (t, x x ) (3.2) 

ip is the field that annihilates a quark and x the one that creates and antiquark; pi = — iV Xi and 
Si — (Ti/2. Another option to express the degrees of freedom is to represent the quark-antiquark 
pair by a wavefunction field 'J (that is to project the theory to the one heavy quark-one heavy 
antiquark sector) 

*(t,x 1( x 2 ) Q/ 3 ~ ^ a (t,xi)xJ(t,x 2 ) ~ -^<5 Q/3 Vv(t,xi)xU*> x 2) + ^T^T a pa ^{t, ^x)x% x 2 ) (3.3) 

Now the Lagrangian has the form (mi is the mass of the heavy quark and m 2 the mass of the heavy 
antiquark, later on we will mainly focus in the equal mass case mi = m 2 = m) 

Lnrqcd = J d 3 xid 3 x 2 Tr |*t( ijXljXa ) f Wo + EL + + . . )j *(t, Xl , Xa )| - 

-J d 3 x -G; u (x)G^ a (x)+ J d 3 x Y,Ux)iJpqi{x)+--- (3.4) 



L pot = J d 3 xid 3 x 2 Tr xi, x 2 )y(r,pi,p2, Si, S 2 )(us gluon fields) xi, x 2 )} (3.5) 

where the dots represent higher order terms in the 1/m expansion and 

iD *(t,xi,x 2 ) = i9 *(i,xi,x 2 ) - 5 A (t,Xi)*(t,xi,x 2 ) + *(t,xi,x 2 ) 3 ylo(t,x 2 ). (3.6) 

The gluon fields can be enforced to be ultrasoft by multipole expanding them in the relative coordi- 
nate r (we define the center of mass coordinates by R = (xi +x 2 )/2 and r = xi — x 2 ), the problem is 
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that the multipolc expansion spoils the manifest gauge invariance of the Lagrangian. The gauge in- 
variance can be restored by decomposing the wavefunction field into (singlet and octet) components 
which have homogeneous gauge transformations with respect to the center of mass coordinate 

*(t,x 1 ,x 2 ) = P[e i9 ^ 1A ' dx ] S(r,R,t) + P[e^.^ lA - dx ] 0(r, R, t) P[e ig & A ' dx ] = 

= Up (xi, R) (S(r, R, t) + 0(r, R, t)) U P (R, x 2 ) (3.7) 

with 

Up (x 1; R) = P[e ig ^ A (*> x )' dx ] (3.8) 
and the following color normalization for the singlet and octet fields 

S = 51 c /v/iVc O = O a T7V7> (3.9) 

Arranging things that way, the lagrangian density (at order p 3 /m 2 ) reads 

CpNRQCD = j d 3 r Tr | S+ (id - h s (v, p, P R , Si, S 2 )) S + O f (<D - h (r, p, P R , S x , S 2 )) o| + 
+V A (r)Tr {O+r • 5 E S + S f r • 5 E 0} + Ysjp. T r {O f r • , 9 E O + O f Or • , 9 E} - 

n f 

-k^'^ + ^ft*, (3.10) 



i=l 



where 



h s (v, p, P R , Si, S 2 ) = + + y s ( r; P; p R; Sl , S 2 ) (3.11) 

TOrod ^ mtot 

h„(r, p, P R , Si, S 2 ) = + + y o(l . ; p , p R; Sl , S 2 ) (3.12) 

m rod 2 m to t 



and 



D O = id O - g[A (R, t), O] Pr = -iD R m rod = miT " 2 m to t = mi + m 2 (3.13) 

m to t 

E 4 = G z0 and B' = — ey^G- 3 /2 are the chromoelectric and chromomagnetic fields, respectively. 
Some of the Feynman rules arising from this Lagrangian are displayed in appendix^] When written 
in terms of these singlet and octet fields, the power counting of the pNRQCD Lagrangian is easy to 
establish. Since the Lagrangian is bilinear in these fields we have just to set the size of the terms 
multiplying those bilinears. The derivatives with respect to the relative coordinate and 1/r factors 
must be counted as the soft scale and the time derivatives, center of mass derivatives and fields for 
the light degrees of freedom must be counted as the ultrasoft scale. The a s that come from the 
matching from NRQCD must be understood as a s (l/r) and the ones associated with light degrees 
of freedom must be understood as a s (E). 

It is not that one form of the Lagrangian is preferred among the others, but the different forms of 
writing the Lagrangian are convenient for different purposes. In principle it is possible to go from one 
form of the Lagrangian to the others; as an easy example consider the leading order Lagrangian (in 
a s and in the multipole expansion) in the static limit (m — > oo) written in terms of the wavefunction 
field 

L pN RQCD = J d 3 X!d 3 x 2 Tr {¥t(i jXljXa )(jD )$(t )XljXa )} + 
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+ / d 3 Xl d 3 x 2 — ^ f Tr(r *t( t)Xl)X2 ) T «^,(t )Xl)X2 ))_ 

J |xi - x 2 | 

- J d 3 x -G° ut (x)G'»" , (x) + I d 3 x J2Qi^)iPli(x) (3.14) 

i— 1 

we will forget about the last line in the equation above, since it remains the same. Now we introduce 
the singlet and octet fields, and take into account that at leading order in the multipole expansion 
the Wilson lines arc equal to one, to obtain 

J d 3 Kd 3 r Tr { (S f + O f ) iD (S + O)} + J d 3 Rd 3 r ^Tr {T a (S f + O f ) T a (S + O)} (3.15) 

now, since iDo(S + O) = ido(S + O) — g [Ao, O] and taking into account that the trace of a single 
color matrix is zero, we obtain from the first term in i|3.15|) 



Tr {S+iflfaS + OUDoO} (3.16) 

and from the second term 

— Tr {T a S^T a S + T a O f T a O} = — Tr Icf&S - ^C^ol (3.17) 

which gives us the static pNRQCD Lagrangian at leading order written in terms of singlet and octet 
fields 

L pNRQ cd = J d 3 Rd 3 r Tr js* (ido + S + 0+ (iD oX 

- f d 3 R^G a ^G^ a + f d 3 Kj2qiip qi (3.18) 

While this procedure is relatively simple at leading order, in general it is more convenient to construct 
each form of the pNRQCD Lagrangian independently (by using the appropriate symmetry arguments 
and matching to NRQCD). 

Note that, as mentioned before, the usual quantum mechanical potentials appear as matching 
coefficients of the effective theory. Renormalization group improved expressions for the potentials 
can then be obtained piTU). 



3.1.2 Strong coupling regime 

In this situation (where, remember, |p| > Aqcd ^> E) the physics at the scale of the binding energy 
(which is in what we are interested) is below the scale Aqcd- This implies that QCD is strongly 
coupled, which in turn indicates that is better to formulate the theory in terms of hadronic degrees 
of freedom. Hence we have, unavoidable, to rely on some general considerations and indications from 
the lattice data to identify the relevant degrees of freedom. Therefore we assume that a singlet field 
describing the heavy quarkonium state together with Goldstone boson fields, which are ultrasoft 
degrees of freedom, are the relevant degrees of freedom for this theory. For this assumption to hold, 
we have to consider that there is an energy gap of order Aqcd from the ground state energy to 
the higher hybrid excitations (that is states with excitations of the gluonic spin), which seems to be 
supported by lattice data, and also that we are away from the energy threshold for the creation of 
a heavy- light meson pair (in order to avoid mixing effects with these states). If one forgets about 
the Goldstone boson fields (switch off light fermions), as it is usually done, we are left with just 
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the singlet field and the theory takes the form of the potential models. In that case the pNRQCD 
Lagrangian is given by 

£ P NRQCD= J d 3 R J d 3 rS^(id -h s (xi,x 2 ,pi,p 2 ,Si,S 2 ))S (3.19) 

with 

2 2 

Mxi,X2,pi,P3,Sl,S 2 ) = 7r 3 " + ^ 2 - + V S (xi,X2,pl,P2,Si,S 2 ) (3.20) 

ZJJli ZTTL2 

The potential V s is now a non-perturbative quantity (the different parts of which can be organized 
according to their scaling in m). The matching procedure will give us expressions for the different 
parts of the potential in terms of Wilson loop amplitudes (which in principle could be calculated on 
the lattice or with some vacuum model of QCD). When considering annihilation processes (in which 
case, obviously, mi = m 2 = m), these expressions translate into formulas for the NRQCD matrix 
elements. Hence, in the strong coupling regime, the NRQCD matrix elements can be expressed 
in terms of wave functions at the origin and a few universal (that is bound state independent) 
parameters. A list of some of the pNRQCD expressions for the matrix elements can be found in 
appendix [Cl 

In the process of integrating out the degrees of freedom, from the scale m to the ultrasoft scale, 
new momentum regions may appear (which were not present in the weak coupling regime, since 
now we are also integrating Aqcd)- It turns out that the intermediate three momentum scale 
y/mAQCD is also relevant (it give contributions to loop diagrams where gluons of energy Aqc d are 
involved. Note that WttiAqcd is the three momentum scale that corresponds to the energy scale 
Aqcd)- Hence, effects coming from this intermediate scale have also to be taken into account for 
the matching in the strong coupling regime |llj . 

To establish the power counting of this Lagrangian we have to assign the soft scale to derivatives 
with respect to the relative coordinate and 1/m factors, and the ultrasoft scale E to time derivatives 
and the static Vs . By definition of strong coupling regime a s evaluated at the scale E must be taken 
as order one. If we want to stay in the most conservative situation we should assume Aqcd ~ mv, 
in which case a s {\/r) ~ 1. Expectation values of fields for the light degrees of freedom should be 
counted as Aqcd to the power of their dimension. 



3.2 Soft-Collinear Effective Theory 

The aim of this theory is to describe processes in which very energetic (collinear) modes interact with 
soft degrees of freedom. Soft-Collinear Effective Theory (SCET) can thus be applied to a wide range 
of processes, in which this kinematic situation is present. Those include exclusive and semi-inclusive 
B meson decays, deep inelastic scattering and Drell-Yan processes near the end-point, exclusive and 
semi-inclusive quarkonium decays and many others. 

Generally speaking, any process that contains highly energetic hadrons (that is hadrons with 
energy much larger than its mass), together with a source for them, will contain particles (referred 
to as collinear) which move close to a light cone direction n/\ Since these particles are constrained 
to have large energy E and small invariant mass, the size of the different components (in light 
cone coordinates, — (hp)n^ /2 + pj_ + (np)n fl /2) of their momentum p is very different; typically 
hp r~j E, p± r~j EX and np ~ EX 2 , with A a small parameter. It is of this hierarchy, hp ^> p± ^> np, 
that the effective theory takes advantage. Due to the peculiar nature of the light cone interactions, 
the resulting theory turns out to be non-local in one of the light cone directions (as it has been 
mentioned in the introduction of the thesis). 
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Unfortunately there is not a standard notation for the theory. Apart from differences in the 
naming of the distinct modes, there are basically two different formalisms (or notations). The 
one originally used in |ol 112). which uses the label operators (sometimes referred to as the hybrid 
momentum-position space representation) and the one first employed in [7|, which uses light-front 
multipole expansions to ensure a well defined power counting (this is sometimes referred to as 
the position space representation) 2 . The two formalisms are supposed to be completely equivalent 
(although precise comparisons are, many times, difficult). 

The modes one need to include in the effective theory depend on whether one want to study 
inclusive or exclusive processes. The resulting theories are usually called SCETi and SCETn, re- 
spectively. When one is studying an inclusive process, collinear degrees of freedom with a typical 
offshellness of order ^JEKqcd are needed. While in an exclusive process the collinear degrees of 
freedom in the final state have typical offshellness of order Aqcd] the simultaneous presence of two 
type of collinear modes must then be taken into account in the matching procedure from QCD to 
SCETn. We will briefly describe these two theories in turn, in the following subsections. In this 
thesis we will be mainly using the SCETi framework (consequently the peculiarities and subtleties 
of SCETn will just be very briefly mentioned). 

3.2.1 SCETi 

This is the theory containing collinear (p^ 1 — (fip,p±, np) ~ (I, A, A 2 )) and ultrasoft (p^ ~ (A 2 , A 2 , A 2 )) 
modes 3 (for some applications collinear fields in more than one direction could be needed), where 
the final collinear states have virtualities of order EAqcd- The theory was first written in the 
(sometimes called) label or hybrid formalism Within that approach the large component of 

the momentum p is extracted from the fields (it becomes a label for them) according to 

<}>(x)=J2e- ii>x <f>n, P (3.21) 

where p = p + k and p contains the large components of the momentum. In that way rip and p± 
have become labels for the field. Derivatives acting on (j) ntP will just give contributions of order A 2 . 
Then the so called label operators V are introduced. Those operators, when acting on the effective 
theory fields, give the sum of large labels in the fields minus the sum of large labels in the conjugate 
fields. We have, therefore 

/ (V) « gi • • • <f> n , Pl ...)=/ (n Pl + ■.. - n qi ■ ■ ■) • • • d> n , Pl • • •) (3.22) 

an analogous operator is defined for the transverse label With that technology, building blocks 
to form invariant operators (under collinear and ultrasoft gauge transformations) can be constructed. 
A scaling in A is assigned to the fields in the effective theory, such that the action for the kinetic 
terms counts as A . The scaling for the various fields is summarized in table l3~2l The leading 

2 Note that the multipole expansions used in the position space representation are, to some extent, similar to the 
ones used in pNRQCD, while the hybrid representation is, not surprisingly, closer to the so called vNRQCD formalism. 
In any case the main difference between the vNRQCD and pNRQCD approaches is the way the soft and ultrasoft 
effects are disentangled. While vNRQCD introduces separate fields for the soft and ultrasoft degrees of freedom at 
the NRQCD level, the pNRQCD approach integrates out the soft scale producing thus a chain of effective theories 
QCD— >NRQCD— >pNRQCD, so that the final theory just contains the relevant degrees of freedom to study physics 
at the scale of the binding energy. In that sense any version of SCET is closer to vNRQCD than to pNRQCD, since 
separate (and overlapping) fields are introduced for the soft and collinear degrees of freedom (which probably is more 
adequate in this case) . 

3 Be aware that the terminology for the different modes varies a lot in the literature. One should check the 
terminology used in each case to avoid unnecessary confusions (this is also true for SCETn). 
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Fields 


Scaling 


collinear quark £ 


A 


ultrasoft quark q 


A 3 


ultrasoft gluon A^ s 


A 2 


collinear gluon (nA niq , A^ q ,nA n ^ q ) 


(1,A,A 2 ) 



Table 3.2: A scaling of the fields in SCET^ 



order Lagrangian for the SCET is then derived. This leading order (in the power counting in A) 
Lagrangian is given by 



£ c = W jmD + gnA n . q + (f ± + g4l q ) W±W^ (f ± + g^ q ,) \ ^ n>p 



(3.23) 



in that equation £ are the fields for the collinear quarks, A arc the gluon fields, the covariant 
derivative D contains ultrasoft gluon fields and W are collinear Wilson lines given by 



W 



perm. 



(3.24) 



where the label operator acts only inside the square brackets. We can see that couplings to an 
arbitrary number of nA n _ q gluons are present at leading order in A. The Feynman rules arising from 
this Lagrangian are given in appendix iBl 

Subsequently power suppressed (in A) corrections to that Lagrangian were derived. This was 
first done in 0Q3], where the position space formalism for SCET was introduced. In the position 
space formalism, the different modes present in the theory are also defined by the scaling properties 
of their momentum. But the strategy to construct the theory consists now of three steps. First 
one performs a field redefinition on the QCD fields, to introduce the fields with the desired scaling 
properties. Then the resulting Lagrangian is expanded, in order to achieve an homogeneous scaling 
in A of all the terms in it. This step involves multipole expanding the ultrasoft fields in one light 
cone direction, according to 

(f) us (x) = (j) us (x-)+[x±d(f> us ] (x-)+-nx [nd(j) us \ i x -) + ^ [x^x^d^d 11 '(p us ] (x_)+C (A 3 (/>„.,) (3.25) 

where x_ = l/2(nx)n. And finally the last step consists in a further field redefinition which restores 
the explicit (collinear and ultrasoft) gauge invariance of the Lagrangian (which was lost by the 
multipole expansions). With that procedure the Lagrangian for SCET up to corrections of order A 2 
(with respect to the leading term H3.23fl ) was obtained. Later on this power suppressed terms were 
also derived in the label formalism |14|. 

Note that the purely collinear part of the Lagrangian is equivalent to full QCD (in a particular 
reference frame). The notion of collinear particle acquires a useful meaning when, in a particular 
reference frame, we have a source that create such particles. 
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3.2.2 SCETn 

This is the theory that describe processes in which the collinear particles in the final state have 
virtualities of order A.q CD . The simultaneous presence of two kinds of collinear modes must be taken 
into account in this case. We will have hard-collinear modes, with a typical scaling p M ~ (1, A, A 2 ) 
and virtuality of order EKqc d (these correspond to the collinear modes of the previous section, in 
SCETi) and collinear modes, with a typical scaling p^ ~ (1, A 2 , A 4 ) and virtuality of order Aq C£) ; 
together with ultrasoft modes with scaling p^ ~ (A 2 , A 2 , A 2 ). 

In the final effective theory (SCETn) only modes with virtuality 0(Aq CD ) must be present. The 
contributions from the intermediate hard-collinear scale must then be integrated out in this case. 
This can be done with a two step process, where first the hard scale E is integrated and one ends up 
with SCETj. Then the hard-collinear modes are integrated and one is left with an effective theory 
containing only modes with virtuality of order Aq CD . SCETn is therefore much more complex than 
SCETi. In particular one of the most controversial issues is how one should deal with end-point 
singularities that may appear in convolutions for the soft-collinear factorization formulas. Those can 
be treated, or regulated, in several different ways. If one works in dimensional regularization in the 
limit of vanishing quark masses a new mode, called soft-collinear messenger |15j . must be introduced 
in the theory. It provides a systematic way to discuss factorization and end-point singularities. 
Alternative regulators avoid the introduction of such a mode. Although this is clear now, to what 
extent the messenger should be considered as fundamental in the definition of the effective theory 
or not is still under debate. 



Chapter 4 

The singlet static QCD potential 



In this chapter we will calculate the logarithmic fourth order perturbative correction to the static 
quark-antiquark potential for a color singlet state (that is the sub-leading infrared dependence). 
This work appears here for the first time. It will later be reported in |16|. 



4.1 Introduction 

The static potential between a quark and an antiquark is a key object for understanding the dynamics 
of QCD. The first thing almost every student learns about QCD is that a linear growing potential 
at long distances is a signal for confinement. Apart from that, it is also a basic ingredient of the 
Schrodinger-like formulation of heavy quarkonium. What is more, precise lattice data for the short 
distance part of the potential is nowadays available, allowing for a comparison between lattice and 
perturbation theory. Therefore the static potential is an ideal place to study the interplay of the 
perturbative and the non-perturbative aspects of QCD. 

The quark-antiquark system can be in a color singlet or in a color octet configuration. Which 
will give rise to the singlet and octet potentials, respectively. Both of them are relevant for the 
modern effective field theory calculations in the heavy quarkonium system. Here we will focus in 
the singlet potential. 

The perturbative expansion of the singlet static potential (in position space) reads 
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the non-logarithmic third order correction 03 is still unknown. The form of the logarithmic term 
in (|4.1I) corresponds to using dimensional regularization for the ultrasoft loop (which is the natural 
scheme when calculating from pNRQCD). 1 

We will calculate here the logarithmic fourth order correction to the potential. Since this calcu- 
lation follow the same lines as that of the third order logarithmic terms, we will briefly review it in 
the next section. 

4.2 Review of the third order logarithmic correction 

The leading infrared (IR) logarithmic dependence of the singlet static potential was obtained in 
jS] by matching NRQCD to pNRQCD perturbatively. The matching is performed by comparing 
Green functions in NRQCD and pNRQCD (in coordinate space), order by order in 1/m and in the 
multipole expansion. 

To perform that matching, first of all one need to identify interpolating fields in NRQCD with the 
same quantum numbers and transformation properties as the singlet and octet fields in pNRQCD. 
The chosen fields are 

X t( X2 , i)0(x a , x i; t)^( Xl , t) - ^zf\r)S{v, R, t) + ^Z E . s (r)rv ■ 5 E a (R, i)O a (r, R, t) + . . . (4.4) 
for the singlet, and 

X t (x 2 ,t)^(x 2 ,R;0T>(R,x 1 ;i)V'(xi,t) -> \/^ 0) (r)O Q (r, R, t) + 

+ ^Z E ,o{r) r r • 3 E a (R, t)S(r, R, t) + . . . (4.5) 

for the octet, where 

<f>{y, x; t) = P exp jt J ds (y - x) ■ 5 A(x - s(x - y), t) J (4.6) 

The Z in the above expressions are normalization factors. The different combinations of fields in 
the pNRQCD side are organized according to the multipole expansion, just the first term of this 
expansion is needed for our purposes here. Then the matching is done using the Green function 

G = (vac|x t (x 2 )0(x 2 ,a;i)V'(a;i)V' t (2/i)'/ , (yi,y2)x(y2)|vac) (4.7) 

In the NRQCD side we obtain 

Gnrqcd = £ 3 (xi - yi)<5 3 (x 2 - y 2 )(W n ) (4.8) 

where Wu represents the rectangular Wilson loop of figure |4~T1 Explicitly it is given by 

W n = Pcxp^-ig£ dz^A^zn . (4.9) 

The brackets around it in 1)4.8(1 represent an average over gauge fields and light quarks. We are 
interested only in the large T limit of the Wilson loop (to single out the soft scale), therefore we 
define the following expansion for T — > 00 

i HW U ) = u {r) + l ^ + O (Jz) (4-10) 

1 Note that this is not the natural scheme when calculating from NRQCD. In that case one would regulate also the 
potentials in d-dimensions. 
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Figure 4.1: Rectangular Wilson loop. 
(-T/2,r/2) and y 2 = (-T/2, -r/2) 



The corners are xi = (T/2,r/2), x 2 = (T/2,-r/2), yi 





Figure 4.2: Schematic representation of the NRQCD^pNRQCD matching for the static potential 
and the normalization factors, at order r 2 in the multipole expansion. On the pNQRCD side (right) 
the single line represents the singlet field, the double line the octet field, the circled cross the 
0(r) chromoelectric vertex and the thick springy line represents the correlator of chromoelectric 
fields. Remember that those diagrams are representing an expansion in r (and in 1/m) and not a 
perturbative expansion in a s . 



In the pNRQCD side we obtain, at order r 2 in the multipole expansion 2 



Gpnrqcd = ZW(r)<5 d ( Xl - yi)^(x 2 - y 2 )e" 

' 2 dt f dt 1 e^ t - t '^ a) - v ^ ) \v ■ g^%t)^{t,t')v ■ g^it')) 

T/2 J-T/2 



(4.11) 



where the Wilson line 



4>{t,t') = P exp 



ig J dtA (t) 



(4.12) 



which comes from the octet propagator, is evaluated in the adjoint representation. 

Then comparing Gnrqcd with G p nrqcd one obtains the matching conditions for the potential 
and the normalization factors. That matching is schematically represented in figure l4~2l Note that up 
to this point no perturbative expansion in a s have been used yet. We will now evaluate the pNRQCD 
diagram perturbatively in a s . The dependence in a s enters through the Va, V s and V potentials 
and through the field strength correlator of chromoelectric fields. To regulate IR divergences we will 
keep the a s dependence in the exponential on the second line of (|4.11|) uncxpandcd. V a — V s will 
then act as our IR regulator. The tree level expression for Va is given by the NRQCD^pNRQCD 
matching at order r in the multipole expansion. It simply states that Va = 1 at tree level. We then 
need the tree level expression for the correlator of chromoelectric fields. Note that, since afterwards 



2 The superscripts (0) in all those expressions are reminding us that we are in the static limit m — > oo. Since in 
this chapter we are always in the static limit, we will omit them after 14.111 . to simplify the notation. 
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we want to integrate over t and if , we need this gluonic correlator in d dimensions (see the following 
subsections for more details). Inserting the rf-dimensional (d = 4 — 2e) tree level expressions we 
obtain (in the T — > oo limit) 



Gpnrqcd = Z s (r)<5 3 (x! - yi)<S 3 (x 2 - y 2 ) 



~2N~ r — —7^ 



,-irv s (r). 



^c-lOsM 2 



27V r 




r 2 (V - V s y - - log -2 + const. (4.13) 



where we can explicitly see that V a — V s acts as our IR regulator in the logarithms. The last line 
in (|4.13|) only affects the normalization factor (and not the potential) and will be omitted in the 
following. The V — V s that appears outside the logarithms must now be expanded in a s (since are 
no longing acting as IR regulators). Therefore wc have obtained the result 



GpNRQCD = Z s (r)5 3 (xi - yi)<5 3 (x 2 - y 2 )e 



-iTV e (r). 



Note that the a s coming from the potentials are evaluated at the soft scale 1/r, while the a s 
coming from the ultrasoft couplings is evaluated at the scale /i 3 . The ultraviolet divergences in 
that expression can be re-absorbed by a renormalization of the potential. Comparing now the 
perturbative evaluations of Gnrqcd and G p nrqcd (obviously the same renormalization scheme one 
has used in the evaluation of the pNRQCD diagram must be used in the calculation of the Wilson 
loop in NRQCD) we obtain 

V s (r,n) = KM)two- 1 oo P s - ^f^^f y ^ln(ru) (4.15) 

when (uo(?"))two-loops is substituted by its known value we obtain the result i|4.1f> . quoted in the 
introduction of this chapter. The calculation explained in this section has thus provided us the 
leading IR logarithmic dependence of the static potential. In [S] the cancellation of the IR cut-off 
dependence between the NRQCD and pNRQCD expressions was checked explicitly, by calculating 
the relevant graphs for the Wilson loop. That was an explicit check that the effective theory 
(pNRQCD) is correctly reproducing the IR. 

In the following section we will use the same procedure employed here to obtain the next-to- 
leading IR logarithmic dependence of the static potential. That is the logarithmic contribution 
to the potential (which is part of the N 4 LO, in a s , correction to the potential). 



4.3 Fourth order logarithmic correction 

In the preceding section no perturbative expansion in a s was used until the paragraph following 
equation (|4.11() . Therefore (|4.11|) is still valid for us here and will be our starting point (note that 
contributions from higher order operators in the multipole expansion have a suppression of order a 2 
with respect to the second term of (|4.11ll . and therefore are unimportant for us here). We need to 

3 Obviously the distinction is only relevant at the next order, that is next section. 
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Figure 4.3: NRQCD diagrams for the leading order matching of Va- The solid lines represent the 
quark and the antiquark, the dashed line represents an A gluon. 




Figure 4.4: Sample NRQCD diagram for the next-to- leading order matching of Va- 

calculate the a s correction to the evaluation of the diagram in the preceding section. Remember that 
the dependence in a s enters through Va, V s , V and the correlator of gluonic fields, therefore we need 
the 0(a s ) correction to all this quantities. That terms will be discussed in the following subsections 
in turn. Then in subsection 14.3.31 we will obtain the fourth order logarithmic correction to the 
potential. Again we will regulate IR divergence by keeping the exponential of V — V s unexpanded. 

4.3.1 0(a s ) correction of Va, V s and V Q 

The 0(a s ) corrections to V s and V are well known. They are given by 

V s = -^a s (l/r)(l + (a 1 +2 lE p )^^\ (4.16) 

V = (^-C f y-a s (l/r)(l + (a 1 + 2 lE (3 )^^j (4.17) 

as it has already been reported in the introduction of this chapter. The mixed potential Va can be 
obtained by matching NRQCD to pNRQCD at order r in the multipole expansion (we are thinking 
in performing this matching for Va in pure dimensional rcgularization). At leading order in a s we 
have to calculate the diagrams shown in figure fOl They give the tree level result Va = 1. The first 
corrections to this result are given by diagrams like that of figure |4~4*I We can clearly see that the 
first corrections are 0{a 2 s ) 

V A = l + 0(a 2 s ) (4.18) 

and therefore unimportant for us here. 

4.3.2 0(a s ) correction of the field strength correlator 

The 0(a s ) correction to the QCD field strength correlator was calculated in |5T]. Let us review here 
that result and explain how we need to use it. 
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Figure 4.5: Leading order contribution to the field strength correlator. The gluonic string is repre- 
sented by a double line. 

The two-point field strength correlator 

V^{z) = (vac| T {G^(y)Ve^ abCzT K d ° A ° ("+") G b Xu (x) } |vac) (4.19) 

can be parametrised in terms of two scalar functions T)(z 2 ) and T>i(z 2 ) according to 
T> f iu\u{z) = (g^xgvuj - g^gvx) ip{z 2 ) + T>i(z 2 )) + 

dT>i(z 2 ) 

Zuj + gvujZfj,z\) — ^2 — (4.20) 

where z = y — x. In (|4.11() x and y just differ in the time component, so z = t — t 1 for us. Furthermore, 
we are interested in the chromoelectric components, so we need the contraction 

V mo (z) = -(d - 1) (v(z 2 ) + V 1 (z 2 ) + ^qP ^ (4-21) 

The tree level contribution is given by the diagram shown in figure l4~5*f . the result is 

Vf{z 2 ) = tf\Nl l) J[ 2 e ~_t V^{z 2 ) = (4.22) 

The next-to-leading (0(a s )) contribution is given by the diagrams in figure 1431 Here we need the 
expression in d dimensions. The d-dimensional result for the a s correction is |22| 

V^\z 2 ) = ^ 2 c -X)^^Y\l-e)(^^ g(e) (4.23) 

2><V) = N^-D^^T^l-e)^^^) (4.24) 

with 

= 2.3 + 2(1-^(2.-1,2.-2)^-6^ + 8.-3 
yy ' e(2e 2 - 5. + 3) y ' 

-6. 3 + 17. 2 -18. + 6 2(l-. + . 2 )i?(2.-l,2.-2) + fcW 

5l(£) = . 2 (2. 2 -5. + 3) + 7(2^3) (426) 

Since the external points, x and y, are fixed in this calculation, the divergences we will encounter in 
D mo (coming from the expressions above) should be canceled by the vertex and (gluon and octet 

4 Note a slight change in notation in the diagrams with respect to 1211 . We represent the gluonic string by a double 
(not dashed) line and we always display it (also when it reduces to S ab ). 

5 We are indebted to Matthias Jamin for sharing the d-dimensional results for the field strength correlator with us. 
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Figure 4.6: Next-to-leading order contributions to the field strength correlator. The gluonic string is 
represented by a double line. The shaded blob represents the insertion of the 1-loop gluon self-energy. 
Symmetric graphs are understood for (c) and (d). 
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Figure 4.7: 0(a s ) covmterterm diagrams for the chromoelectric correlator. The gluonic string (which 
comes from the octet propagator) is represented by a double line. 

field) propagator counterterms. The counterterm for the vertex is zero, since as we have seen in 
the previous subsection the first correction to Va is of order oq. The counterterm for the gluon 
propagator is the usual one in QCD. The counterterm for the octet propagator coincides with the 
counterterm for the quark propagator in Heavy Quark Effective Theory but with the quark in the 
adjoint representation. We can represent the counterterm contributions by the diagrams of figure 
14.71 We have checked that when we compute 2?ioio the divergence coming from the first diagram in 
figure EHJI is canceled by the counterterm of the gluon propagator. That diagram l4.6b does not give a 
divergent contribution, as it should. And that when we add the remaining diagrams the divergence 
we obtain is exactly canceled by the counterterm of the octet propagator. 
The contributions of the counterterms are given by 



where the first 1/e in the square bracket corresponds to the octet propagator and the second one 
to the gluon propagator. Then the total d-dimensional result (including the contributions from the 
counterterms) for the a s correction to the chromoelectric correlator is 



£>£] = -(d- 1) (v {1 Hz 2 ) + (-1 + 2e)v[ 1) (z 2 ) +V c - t -(z 2 ) + (-1 + e)^^ 2 )) (4.29) 



which no longer have 1/e poles. 

4.3.3 Calculation of the fourth order logarithmic correction 

The results of the two preceding subsections provides us all the necessary ingredients to compute 
the next-to-leading IR dependence of the static potential. We have to evaluate the second term in 
the parenthesis of (f 4 . 1 1|> at next-to-leading order. Let us define 



= -if Vl(r) / dt / dt> e -«('-'')(Vo-v.) (r . 5 E a (t)C' j (M> ■ flE 6 (i')) (4.30) 



First we will consider the contribution we obtain when we insert the a s correction l|4.29|l to the field 
strength correlator. We have just to perform the integrations over t and t' . To do that we change 
the integration variables to t + t' and t — t' = t-. The integral over the sum just gives us a factor 
T — t— . The T term will give a contribution to the potential and the i_ term a contribution to 
the normalization factor (that is unimportant for us here, it will be omitted in the following). The 
remaining integral over t- can be done by using 




(4.27) 



(4.28) 




■T/2 




(4.31) 
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The result we obtain is then (in the T — > oo limit) 



0&>l =-iT(^-) aKl/r)^!^ 1 -. 
*< EE >\o( as ) \ 7T / sW ' 2 8 r 

Ci log h 6 2 log h const. (4.32) 



e /i /i 



with 



A = — 



1 (2n s 11 



24 V 3iV c 3 
6\ 3iV c 3 

C 2 = 12^-149 j (4.33) 

We get another contribution to ( r ' when we use the leading order expression for the chromoelectric 
correlator (then we arrive at (|4.13l) ) and then insert the next-to-leading order correction to V a — V s . 
This contribution is given by 

Mr 2 ) a s (^)a A J\lr)Nl-lC\ 1 . „ n .fl , fV -V s \ 2 \ 

(4.34) 

The ultraviolet divergences we encounter in expressions (|4.32|l and (|4.34|) can again be re-absorbed 
by a renormalization of the potential. Finally we get another contribution that comes from changing 
a s (u) to a s (l/r) in equation l|4.14|l . after renormalization (we want all the a s evaluated at the scale 
1/r in the potential) . It is given by 

G^\ = -ySZWm* * ~ ^logMlog (YozVL) (4.35) 
^ 1 / r 48tt 2 2N C \ fi J ■ ' 

We see that the log 2 ((V r Q — V s )/fi) and the log(r/i) log((V r c — V s )/(j,) terms appear with the right 
coefficients to form, together with the double IR logarithms that would come from the NRQCD 
calculation of the Wilson loop, an IR cut-off independent quantity for the matching coefficient (as it 
should). Moreover the coefficient for the double logarithm we have obtained here (which remember 
came from the correction to the gluonic correlator) coincides with what one obtains expanding the 
renormalization group improved static potential of • These two facts are checks of our calculation. 

We have therefore obtained the a|? log rzz (and a|? log 2 r/x) terms of the singlet static potential 6 . 

V s (Q) (r) = (EqETJ- 

(..smy l|! c s + ?„,) ^ T „ _ (4 . 36) 



4?r 



16^ 2 Ci ^ai + 2 1E (3 - ^ (20n y - C a (12tt 2 + 149))) logr/i (4.37) 



6 Note added: It is understood that (when renormalizing the potential) we have used the scheme where 1/e — "fE + 
log n is subtracted. And this has been implemented by redefining fj, 2 — > fi 2 e lE /n where applicable. Therefore one 
of the a B in the third order correction in 14.11 is understood to be in this scheme, whereas the remaining ct B are 
understood to be in the MS scheme. Also we have chosen the scheme where only the log rfi terms that compensate 
an infrared log((V — V s )/fJ,) are displayed in the potential. 
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4.4 Discussion 

In view of the possible future construction of an e + — e~ linear collider (the aim of which will be 
the study of the possible new particles LHC will discover), much theoretical effort is being put in 
the calculation of t — i production near threshold. The complete second order (N 2 LO) corrections 
are already computed. The second order renormalization group improved expressions (N 2 LL) are 
under study (several contributions are already known) (23 ISI I2H HOI CHI • Given the extremely good 
precision that such a new machine could achieve, the third order corrections are also needed. These 
third order corrections (N 3 LO terms) are being computed at present, by several different people. 
This is a gigantic project that requires the use of state-of-the-art calculational and computational 
techniques [2E1 EH E0I OH] • Once those third order corrections are completed, the corresponding- 
third order renormalization group improved expressions (N 3 LL) will also be needed (to achieve the 
desired theoretical precision in the calculation). Just let us mention that the results presented in 
this chapter will be a piece of these N 3 LL computations. 
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Two loop SCET heavy-to-light 
current anomalous dimension: rtj 
terms 



In this chapter we will calculate the two loop n / terms of the anomalous dimension of the leading 
order (in A) heavy-to- light current in Soft-Collinear Effective Theory (SCET). The work presented 
in this chapter, although mentioned in appears here for the first time. The calculation of 

the complete two loop anomalous dimension will appear elsewhere |34| . 

5.1 Introduction 

The heavy-to-light hadronic currents Jh ac i = q^b (b represents a heavy quark and q a light quark) 
appearing in operators of the weak theory at a scale ~ mj can be matched into SCETj . The 
lowest order SCET hadronic current is not J^ d ET = C(fj,)^Th, but rather 

^had ET = c o (np, (i) £n,z>r/i + ci (np, nqi,(Ji) £,n, P (gnA n>qi )Th H (5.1) 

That is, an arbitrary number of nA nA gluons can be added without suppression in the power 
counting. Here £ and A are the fields for the collinear quarks and gluons in the effective theory, 
respectively; h is the field for the heavy quark in HQET. Collinear gauge invariance relates all these 
operators and organize the current into the (collinear gauge invariant) form 

4£i ET = Ci{v,nP)xn,pTh (5.2) 

where 

X = (5.3) 

and W is a collinear Wilson line (see section I3~2|) . We can then run the Wilson coefficients down in 
SCET. Note that it is enough to consider the simpler current £r/i, because collinear gauge invariance 
relates them all. This was done at one loop in 6 . The result obtained there was 1 

1 The coefficients for the different Dirac structures mix into themselves. There is no operator mixing at this order. 
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Figure 5.1: Feynman rule for the O(A ) SCET heavy-to-light current. The double line is the heavy 
quark. The dashed line is the collinear quark. Springy lines with a line inside are collincar gluons. 

i 1 + f 



vk + irj 2 



= igvH a 



Figure 5.2: HQET Feynman rules. The heavy quark is represented by a double line. The springy 
line represents the gluon. 
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(5.5) 



Z is the current counterterm in the effective theory, 7 is the anomalous dimension (P is the total 
outgoing jet momentum). Here we will calculate the 2 loop rif corrections to this result. 



5.2 Calculation of the rif terms 

The effective theory diagrams that are needed to compute the n f terms of the two loop anomalous 
dimension are depicted in figure 15.31 We will perform the calculation in d = 4 — 2e dimensions. 
To distinguish infrared (IR) divergences from the ultraviolet (UV) ones, we will take the collinear 
quark off-shell by setting p± = and the heavy quark with residual momentum oj. This will regulate 
the IR divergences of all the diagrams. We will work in Feynman gauge. The gluon self-energy is 
the same as in QCD (for both the collinear and the ultrasoft gluons), it is given in figure IOI The 
Feynman rules which are needed to compute the diagrams are given in figure 15.11 (for the current) , 
figure l5~2l ( vertex and propagator rules for HQET) and appendix [5] (vertex and propagator rules for 
SCET). For the ultrasoft diagrams the collinear quark propagator simplifies to (s is an ultrasoft 
loop momentum, p is the external collinear quark momentum) 



n(p + s) n(p) 



(p + s) 2 + it] n(p)n(p + s) + irj ns + ^2 + i-q n s + ]~ + irj ns + a + irj 

To further simplify the integrals we will choose oj to be u> — a/2. 

For the evaluation of the ultrasoft graphs we will just need the integrals 

1 1 ( 1 \ P _ 2i t 2 _ /3 r(2-/3-e)r(-2 + 2/3 + 2e) 



(5.6) 



7-{-l) 



(2-n) d ns + a + irj vs + lu + irj \s 2 + ir]J (An) 2 e T(/3) 
■ f 1 dy (2™ + a(l - V)f- 20 - 2 ' V- 2+ ^ ^M)" V ^- 2e r(2 - g - £ )r(-2 + 2(3 + 2e) 

J d y^y + a ^ v)) y - (47r) 2-e m{ _ 1 + (3 + e) 
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1 ^ P 2i(-l)^ (2uj)3 - 2 0-2e r(2 - 13 e)r(-3 + 2(3 + 2e) 



7 (27r) d vs + w + z?7 Vs 2 + Z77 / (47r) 2 - £ r(/3) 

which can be calculated with Feynman and Georgi parameterizations, we have used that 2lu — a in 
the last step of 1)5. 7|) . Using these results we obtain 



(47r) 2 \^np J \ \12t A 36e 2 e \ 27 



1 5 1/8 7T 2 



■^/I-O-O + tI^-t))) ( 5 - 9 ) 



3e d 9e 2 e V 27 2 

-2e 



/p 2 \"„/„/5 1/ 5 25tt' 



(4tt) 2 Va"W V V 3e3 e V 3 36 

M-s+Ks-^))) 

where we have redefined /j, 2 — » /^ 2 e 7E /(47r) (from now on, we will always use this redefinition). 

The evaluation of the collinear graphs requires the integral (which again can be calculated with 
Feynman and Georgi parameterizations) 

d d s n(p—s) 1 / 1 \^ i{-p 2 )~ £ , 2\-/3+i 



(2-7r) rf ns (s -p) 2 + it \s 2 + it) (47r) 2 - e yV 
r(l - e)r(/3 - 1 + e)r(l - p - e) / 1 - e 



r(/3)r(2 - - 2e) V 2 - ^ - 2e / 

plus other integrals which do not involve the n vector and can thus be found in standard QCD books 
(see for instance |35|). Using these results we obtain 



-2c 



Fi ^=(4H^ C \ C -\^^ t K 27 ' 36 
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, 2 16 1 / 176 7T 2 , i > 



-M-*-* + K- j+ £))) 

To compute the anomalous dimension we also need the two loop correction to the collinear 
and heavy quark propagators. Since the correction to the collinear quark propagator just involves 
collinear particles (and not ultrasoft ones), this is the same as in usual QCD. While the correction 
for the heavy quark propagator is that of HQET. The corresponding counterterms are |36l 137) 



z f = i ■ as0/ 



s C f 2 /a s x 2 / / 9 19\ 2 (2 
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Figure 5.3: Effective theory diagrams contributing to the 2 loop n/ terms of the O(A ) heavy-to- 
light current. The double line represents a heavy quark, the dashed line represents a collinear quark 
and the springy lines are gluons (collinear if they have a line inside ultrasoft if not). The circled 
cross is the insertion of the current. The shaded blob represents the one loop insertion of the gluon 
self-energy and the cross the corresponding counterterm. 



<«"«■-«•«•> HI 1 *' -H 

/ 2 (e) - r(e) (B(l - e, 1 - e) + 45(2 - e, 2 - e)+ 
+e(B(l-e,l-e)-4B(2-e,2-e))) 

Figure 5.4: Gluon self-energy graph. The gluons in this graph can be understood either as ultrasoft 
or collinear. 
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With this we obtain the two loop n / part of the counterterm in the effective theory 

z "'°""> - m' \°>^> (h + h + 7M£?)-^\-i-b°z (&)) <«»> 

where P is the total outgoing jet momentum. We can then obtain the two loop n/ terms of the 
anomalous dimension by using the formula 

2 d v 2 (( r a2 s\ dZ , A* 9Z\ 

1= Zdil Z= z{ [- eas - 0o ^)d^ + 2djl) (5 ' 16) 

The result is 

J.2loo r „,, = * T ^I°1 (™ + *L + 5 ,„„ (± 



5.3 Discussion 



There is a lot of theoretical interest in obtaining the O (a 2 s ) corrections to the B — > X s j decay rate. 
To measure this decay rate it is experimentally necessary to put a cut on the energy of the observed 
photon. This cut introduces a new scale in the problem and, consequently induces a possible new 
source of corrections that must be taken into account in the evaluation of the decay rate. The effects 
of this scale can be systematically treated in an effective field theory framework using SCET |32| . 
A factorization formula for the decay rate with the cut in the photon energy which disentangle 
the effects of all these scales, can then be derived in a systematic way. The anomalous dimension 
calculated in this chapter enters in this expression. This formula involves, among many other things, 
a jet function, which describes the physics of the hadronic final state, and a soft function, which 
governs the soft physics inside the B meson. Expressions for the evolution equations of these jet and 
shape functions can be derived. The evolution equation for the jet function involves an anomalous 
dimension j J which can be related to the anomalous dimension for the jet function appearing in 
deep-inelastic scattering (DIS) |221|2H1. The two loop anomalous dimension entering the evolution 
of the shape function was calculated in |39| 2 . The current anomalous dimension that we are dealing 
with in this chapter can be related to the anomalous dimensions for these jet and shape functions 
|3"2"] . Recently the soft and jet functions have been evaluated at two loops This calculation 

confirms the previous results for the two loop anomalous dimension of the soft function (23 EH! 
and provides the first direct calculation (that is, not using the relation with DIS) for the two loop 
anomalous dimension of the jet function. Therefore the last thing that remains to be done is the 
direct two loop calculation of the leading order heavy-to-light current in SCET (in this chapter 
we have presented the rif part of this calculation). This is important since it will ensure that we 
correctly understand the renormalization properties of SCET (at two loops). Given the peculiar 
structure of SCET (much different from other known effective field theories) some subtleties may 
arise here. 



2 The original result in [Ml (from 1992) have been recently corrected in the revised version of the paper (from 
2005). The revised result agrees with an independent calculation 1401 . The discovery of some of the mistakes in 39 
was triggered by the calculations described in this chapter. 



Chapter 6 

Radiative heavy quarkonium 
decays 

In this chapter we study the semi-inclusive radiative decay of heavy quarkonium to light hadrons 
from an effective field theory point of view. As we will see below, the correct treatment of the 
upper end-point region of the spectrum requires the combined use of Non-Relativistic QCD and 
Soft-Collincar Effective Theory. When these two effective theories are consistently combined a very 
good description of the experimental data is achieved. The photon spectrum can then be used to 
uncover some properties of the decaying quarkonia. The contents of this chapter are basically based 
on work published in |42l 1481 14*41 I45| . although some comparisons with new data (not available when 
some of the preceding articles were published) are also presented. 

6.1 Introduction 

Although we will focus on the study of the semi-inclusive radiative decay, the exclusive radiative 
decays have also been addressed in an effective field theory approach. Exclusive decays will be very 
briefly commented in subsection 16.1.21 

6.1.1 Semi-inclusive radiative decays 

Semi-inclusive radiative decays of heavy quarkonium systems (see [2] for a review) to light hadrons 
have been a subject of investigation since the early days of QCD [461 147) . In these references, the 
decay of the heavy quarkonium state to gg^ (and to ggg) is treated in lowest order QCD, in analogy 
with the QED decays of orthopositronium to three photons. This lowest order QCD calculation 
predicted a, basically, linear rise with z (z being the fraction of the maximum energy the photon 
may have) of the photon spectrum. The angular distribution has also been studied in 07j; it 
should be mentioned that this angular distribution is still assumed to be correct, and it is used 
for the comparison of the experimental results with theory and the subsequent extraction of QCD 
parameters 1 . The upper end-point region of the photon spectrum (that was obtained by several 
later experiments @H1 EDI 03) appeared to be poorly described by this linear rise; a much softer 
spectrum, with a peak about z ~ 0.6 — 0.7, was observed instead. A subsequent resummation of 
the leading Sudakov (log(l — z)) logarithms [52], as well as a calculation of the leading relativistic 

1 The recent data in 1481 has allowed for the first time for a check of this prediction for the angular distribution. 
There it is found that data agrees adequately with 1471 
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Figure 6.1: Comparison of the spectrum obtained by the modeling of the octet shape functions 
with data. The dot-dashed line is the (leading order) color singlet contribution alone (color singlet 
model). The dashed line is the direct octet contribution (where the shape functions enter). The solid 
line is the total result. We can clearly see that color singlet model alone is unable to reproduce data 
and that the modeling of the octet shape function has produced a result in complete disagreement 
with data. A is a parameter in the model for the shape functions. Plot from hep-ph/0010217 |6U| . 



corrections |53j (see also |54| 1 . although produced a softening of the spectrum in the upper end-point 
(namely z — > 1) region, were neither able to reproduce the observed spectrum. Instead, the data 
was well described by the model in ref. |55|. where a parton-shower Monte Carlo technique was used 
to incorporate the effects of gluon radiation by the outgoing gluons in the decay. 

This led to some authors to claim that a non-vanishing gluon mass was necessary in order 
to describe the data With the advent of Non-Relativistic QCD (NRQCD) 3 , these decays 
could be analyzed in a framework where short distance effects, at the scale of the heavy quark 
mass m or larger, could be separated in a systematic manner |S7j. These short distance effects are 
calculated perturbatively in a s (m) and encoded in matching coefficients whereas long distance effects 
are parameterized by matrix elements of local NRQCD operators. Even within this framework, a 
finite gluon mass seemed to be necessary to describe data J2HJ. However, about the same time it 
was pointed out that in the upper end-point region the NRQCD factorization approach breaks down 
and shape functions, namely matrix elements of non-local operators, rather than NRQCD matrix 
elements, must be introduced |59|. Early attempts to modeling color octet shape functions produced 
results in complete disagreement with data [HOj (as shown in figure l6~T|l . and hence later authors did 
not include them in their phenomenological analysis. 

Notwithstanding this upper end-point region has received considerable attention lately, as it 
was recognized that Soft-Collinear Effective Theory (SCET) |HH may help in organizing the 
calculation and in performing resummations of large (Sudakov) logs jHSESEJEl]- In fact, the 
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early resummation of Sudakov logarithms |52| has been recently corrected |65| within this framework, 
and statements about the absence of Sudakov suppression in the color singlet channel [HH] have been 
clarified |S3|- These SCET calculations will be explained in the following sections. 

For the T(15) state, the bound state dynamics is amenable of a weak coupling analysis, at least as 
far as the soft scale (mv, v ~ a s (mv) <^ 1, the typical velocity of the heavy quark in the quarkonium 
rest frame) is concerned p7flRj51l(^l7ni[7T1[7!Hl7lfl These calculations can 

most conveniently be done in the framework of potential NRQCD (pNRQCD), a further effective 
theory where the contributions due to the soft and ultrasoft (~ mv 2 ) scales are factorized |80M81l l5] 
(see section ETTJ) . The color octet shape functions can then be calculated combining pNRQCD and 
SCET. This calculation of the octet shape functions in the weak coupling regime will be the subject 
of subsection 16.2.31 

Parallel to all that, shortly after [3], in it was pointed out that a parametrically leading 
contribution had been ignored so far. This was the contribution where the photon is emitted from 
the decay products of the heavy quark (light quarks), and not directly from the heavy quark itself 
(remember that we are always dealing with prompt photons, that is photons that do not come 
from hadronic decays). These type of contributions, called fragmentation contributions, completely 
dominate the spectrum in the lower end-point (namely z — > 0) region. At first it was thought that 
only the gluon to photon fragmentation function appeared in the process; so the radiative decays 
seemed a good place to determine this (yet unknown) gluon to photon fragmentation function; but a 
subsequent investigation [HZ] showed that this was not the case. When considering also the color octet 
contributions, the quark to photon fragmentation function also appeared, and their contributions 
can not be disentangled. 

When all the known contributions to the photon spectrum are taken into account and are con- 
sistently combined, a very good description of the data is now achieved (with no longer need for the 
introduction of a finite gluon mass). This will be explained in detail in section T6. 31 

6.1.2 Exclusive radiative decays 

Exclusive radiative decays of heavy quarkonium have been analyzed in an effective field theory 
framework in jHH]. A combination of NRQCD and SCET is also needed in this case. However, since 
in this case we are dealing with an exclusive process, the existence (and effects) of two different 
collinear scales have to be taken into account. Moreover, the fact that the hadronic final states 
must be composed of collinear fields in a color singlet configuration, causes that only color singlet 
contributions, and not color octet ones, enter in the NRQCD-SCET analysis at leading order (in 
contrast with the situation in the inclusive case, as we will see below). In this case the final result of 
this effective theory analysis agrees with the leading- twist order of previous known results O] |HS] ■ 

We will move now to the study of the semi-inclusive radiative decays, starting in the next section. 

6.2 Effective Field Theory approach to the upper end-point 
region 

The NRQCD framework organizes the radiative decay in the following factorized form 



where H represents a generic heavy quarkonium state and M represents its mass. In that formula 
Ci are the hard matching coefficients, which incorporates short distance effects, and the (H| O t |H) 




(6.1) 
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are the NRQCD matrix elements, which parameterize the long distance effects. 

However, as was already mentioned in the previous section, in the upper end-point region of 
the photon spectrum that standard NRQCD factorization is not applicable This is due to the 
fact that small scales induced by the kinematics enter the problem and have an interplay with the 
bound state dynamics. In order to study this region, one has to take into account collinear degrees 
of freedom in addition to those of NRQCD. This can be done using SCET as it has been described 
in |62l I64| . Using the SCET framework, the decay rate has been expressed in the factorized form 

m 

— = ^2h(M,io,h) / dfc+S'(fc + ,/i)ImJ w (fc + + M(1 - z),n) (6.2) 

where H encodes the short distance effects, J is the so called jet function, which incorporates effects 
at the collinear scale, and S are the ultrasoft shape functions. 

Using the combined counting in SCET (counting in A ~ *J Aqcd/ (2m)) plus NRQCD (counting 
in v), one can see that we have color singlet and color octet operators contributing at the same order. 
More concretely |83l I86j . at 0(A) in the SCET counting we have the 1 5o and 3 Pj octet operators 

E^^U-p^p ( 6 - 3 ) 



E cJ 8,Spj) r^xt P B2A • £a • <tVp (6-4) 

i 

When considering also the v counting (and taking into account the overlap with the 3 Si quarkonium 
state) the two of these operators become 0(u 5 A). Their matching coefficients start at order y/a s {nh)- 
At 0(A 2 ) in the SCET counting and with a matching coefficient starting at order a s (/j,h), we have 
the color singlet operator 

E r W-P A • C l 1,3Sl) Bf l} ( 6 - 5 ) 

i 

when considering also the v counting this operator becomes 0(u 3 A 2 ). Then the octet-to-singlet ratio 
(considering A ~ v) becomes , v ; hence the color octet contributions become as important as 

the color singlet ones if we count a s (fj,h) ~ v 2 ~ 1 — z. 

Before going on, and explaining the calculations in the Effective Field Theory (EFT) approach, let 
us comment on the relation of these EFT calculations with the, phenomenologically very successful, 
model in |55| . Recall that, as mentioned in the introduction of this chapter, in that reference a 
parton-shower Monte Carlo technique was used to incorporate the effects of gluon radiation by the 
outgoing gluons in the decay. When we are at the end-point region of the spectrum, the EFT 
analysis tells us that the decay is organized according to eq. I|6. 2JI . then (as was already explained 
in ref. |64j) the approach in |55) is equivalent to consider that the collinear scale is non-perturbative 
and introduce a model with a gluon mass for the jet function J. When we are away of the upper 
end-point region, the EFT approach tells us that the decay is organized according to Ij6.1|l . Then the 
effects of the gluon radiation modeled in should be incorporated in higher order NRQCD local 
matrix elements (this interpretation is consistent with the fact that the analysis in [55| produced a 
not very large correction to the spectrum for all z, except in the upper end-point region, where the 
effect becomes 0(1)); in any case it is not justified, from an EFT point of view, why one should 
take into account the subset of corrections that are incorporated in |55| and not other ones, which 
in principle could contribute with equal importance. 
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6.2.1 Resummation of the color singlet contributions 

The resummation of the Sudakov logarithms in the color singlet channel has been performed in ref. 

an d in ref. (where the (small) effect of the operator mixing was taken into account). It 
is found that all the logarithms come from collinear physics, that is only collinear gluons appear in 
the diagrams for the running of the singlet operator. 
The resummed rate is given by 
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This result corrects a previous calculation |S2] performed several years ago. 
6.2.2 Resummation of the color octet contributions 

The resummation of the Sudakov logarithms in the color octet channel was performed in ref. |62| . In 
contrast with the color singlet channel, both ultrasoft and collinear gluons contribute to the running 
of the octet operators. The expression for the resummed Wilson coefficients is 



C(x-z) = -^-\6{x-z) 



exp[l gi [g s /3 l/(47r)] + g 2 [g s /3 l/(47r)]] 
r[l - 5l [a s/ V/(47r)] - a s /3 £/(47r) 3 ' 1 [a s/ 3o^/(4^)]] 
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6.2.3 Calculation of the octet shape functions in the weak coupling regime 

In this subsection we will explain in detail the calculation of the octet shape functions. We will start 
by rewriting some of the expressions in the preceding sections in the pNRQCD language, which is 
the convenient language for the subsequent calculation of the shape functions. We begin from the 
formula given in |64| 

^ = z^ImT(z) T{z) =-ij tfxe-*** (V Q (nS)\T{ J^x) J u (0)}\ V Q {nS)) iff (6.15) 

where Jfi{x) is the electromagnetic current for heavy quarks in QCD and we have restricted ourselves 
to 3 Si states. The formula above holds for states fulfilling relativistic normalization. In the case 
that non-relativistic normalization is used, as we shall do below, the right hand side of either the 
first or second formulas in H6.15fl must be multiplied by 2M. At the end-point region the photon 
momentum (in light cone coordinates, q± — q° ± q 3 ) in the rest frame of the heavy quarkonium is 
q = (q+,q-,q±) = (zM, 0,0) with z ~ 1 {My/1 - z < M). This together with the fact that the 
heavy quarkonium is a non-relativistic system fixes the relevant kinematic situation. It is precisely 
in this situation when the standard NRQCD factorization (operator product expansion) breaks 
down |fi9| . The quark (antiquark) momentum in the QQ rest frame can be written as p = (po,p)> 
Po = m + lo,p = 1; Zo?l ^ m - Momentum conservation implies that if a few gluons are produced 
in the short distance annihilation process at least one of them has momentum r — (r+, r_, r±), 
r_ M/2 ; r + ,r± <C M, which we will call collinear. At short distances, the emission of hard 
gluons is penalized by a s {m) and the emission of softer ones by powers of soft scale over M. Hence, 
the leading contribution at short distances consists of the emission of a single collinear gluon. This 
implies that the QQ pair must be in a color octet configuration, which means that the full process 
will have an extra long distance suppression related to the emission of (ultra)soft gluons. The next- 
to-leading contribution at short distances already allows for a singlet QQ configuration. Hence, the 
relative weight of color-singlet and color-octet configurations depends not only on z but also on the 
bound state dynamics, and it is difficult to establish a priori. In order to do so, it is advisable to 
implement the constraints above by introducing suitable EFTs. In a first stage we need NRQCD 
which factors out the scale m in the QQ system, supplemented with collinear gluons, namely 
gluons for which the scale m has been factored out from the components r + , r± (but is still active 
in the component r_). For the purposes of this section it is enough to take for the Lagrangian of 
the collinear gluons the full QCD Lagrangian and enforce r + , r± <C m when necessary. 
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Matching QCD to NRQCD+SCET 

At tree level, the electromagnetic current in (|f>. 1 f>|) can be matched to the following currents in this 

eft jnu 

W = e- i2mX ° (rfeS ^)^) + r & S0)j (8,is 0) (x) + r^j$ pj) (x) + ...)+ h.c. (6.16) 

r (8,^ ) = Jf 8 , 1So) (x) - 

r£uf 5 = ^ + ^ - -^Hp,) (*) = -i^ivv^w (6.i7) 

where n = (n+,ri_,nj_) = (1,0,0) and — e apP o nP ■ These effective currents can be identified 
with the leading order in a s of the currents introduced in |64| (which has already appeared at the 
beginning of this section). We use both Latin (1 to 3) and Greek (0 to 3) indices, is a single 
collincar gluon field here, and eeQ is the charge of the heavy quark. Note, however, that in order to 
arrive at l|6.16|l one need not specify the scaling of collinear fields as M(A 2 , 1, A) but only the cut-offs 
mentioned above, namely r + ,rj_ <C M. Even though the P-wave octet piece appears to be 1/m 
suppressed with respect to the iS-wave octet piece, it will eventually give rise to contributions of the 
same order once the bound state effects are taken into account. This is due to the fact that the 3 5i 
initial state needs a chromomagnetic transition to become an octet 1 S'o, which is a s suppressed with 
respect to the chromoelectric transition required to become an octet 3 Pj. 
T(z) can then be written as 

T(^\ — rriii'aa'PP' . rriS, 1 So)rpaa' , tt(&, 3 Pj) rpaija'i'j' , in -, 8 \ 

\ I ~ U ii' aa 'f3l3' 1 (l3S 1 ) +U aa> 1 (8, 1 S ) + n a ija>i'j> 1 (8,3Pj ) H [0.10) 

where 

rr(l, 3 Sl) _ M ! / r (l, 3 S 1 )p(l, 3 S 1 ) 

11 ii l aa'l3f3' ~ 'l± 1 af3ifi L a'fi'i'v 

rr(8, 1 So) = A"' r (8, 1 S ) r (8, 1 So) 
aa' '/_L ufi 1 a'p 



and 



rr(8, 3 P./) = ^ r (8, 3 -P./) r (8, 3 P./) ,- fil QN 
odja'i'j' 'X afiij cx/ui'j' 

T^'f'(z) = -if ^-^--^(fq^iti^^^^w^o)}!^^ 

T$ So) {z) = -i f d^xe-*-*- 2 ™* (v Q (nS)\T{^ 1So) (x)^ tlSo) (0)}\V Q (nS) 

T$il{{z) = ~i j ' fae-**-*"** (^Q(n5)|T{J^ Pj) (x)V;% ; )(0)}|^ Q (n5)) (6.20) 

In Qfi.lSjl we have not written a crossed term (8, 1 Sq-^Pj) since it eventually vanishes at the order 
we will be calculating. 

Matching NRQCD+SCET to pNRQCD+SCET 



Thanks to the fact that in the end-point region (M ^> My/1 — z ^> M(l — z)) the typical three 
momentum of the heavy quarks is given by 



M . 

— {l-z)-E x ) (<i.2D 
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Figure 6.2: Color octet contributions. • represents the color octet S-wave current, ▲ represents the 
color octet P-wave current. The notation for the other vertices is that of ref. [ST], namely 

-JtB^ {eTl 2° 2) Tr [ ThB ] and ° = 7^^ xTr [ TbE ] ■ Thc solid line represents the singlet field, the 
double line represents the octet field and the gluon with a line inside represents a coUinear gluon. 



we can proceed one step further in the EFT hierarchy. NRQCD still contains quarks and gluons 
with energies ~ ma s , which in the kinematical situation of the end-point (where the typical three 
momentum is always much greater than the typical energy) can be integrated out. This leads to 
potential NRQCD (pNRQCD). 

For the color singlet contributions we have 

(^(n^lTlJ^j^V^^O)}!^^))^ 

— ► 2iV c St, t (x,O,iz;o)St(O,O,0) (vAG\Tr{BfB^}(x)Tr{B^B^ }(0)|VAc\ (6.22) 

The calculation of the vacuum correlator for coUinear gluons above has been carried out in |64| , and 
the final result, which is obtained by sandwiching (|6.22() between the quarkonium states, reduces to 
the one put forward in that reference. 

For the color octet currents, the leading contribution arises from a tree level matching of the 
currents (|6.16(l . 



(8,^0) 



(x) 



V2I>O^(x,0,x )Br(3 



7«J 



v/21> tV'0«?(x, yja: o) 



y=o 



Bl a (x) 



(6.23) 



Sy, Oy and Op are the projection of the singlet and octet wave function fields introduced in |80ll5] to 
their vector and pseudoscalar components, namely S = (Sp + Sya 1 )/^ and O a = {Op + O v l (j l )/\/2. 
Tf = 1/2 and N c = 3 is the number of colors. 



Calculation in pNRQCD+SCET 

We shall then calculate the contributions of the color octet currents in pNRQCD coupled to coUinear 
gluons. They are depicted in figure l6~2l For the contribution of the P-wave current, it is enough 
to have the pNRQCD Lagrangian at leading (non-trivial) order in the multipole expansion given in 
UniEJ. For the contribut ion of the S-wave current, one needs a 1/m chromomagnetic term given in 

m 

Let us consider the contribution of the S-wave color octet current in some detail. We have, from 
the first diagram of fig. 16.21 

32_ 9 /c F x2 , _ f . / \ fd 4 k k 2 



T^s )^) = -ivT\^)Y T F {^Y a M°f J / d 3 x'Co(x')^no(x) /- 



(2tt) 4 k 2 + ie 
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E n -h + ie) x , (M (1 - z) - k + )M -k]_+ie \ -k + E n - h, 



(6.24) 

0,x 



where we have used the Coulomb gauge (both for ultrasoft and collinear gluons). E n < is the 
binding energy [M = 2m + E n ) of the heavy quarkonium, -0no( x ) its wave function, and h the 
color-octet Hamiltonian at leading order, which contains the kinetic term and a repulsive Coulomb 
potential 80, 5 . cp is the hard matching coefficient of the chromomagnetic interaction in NRQCD 
P], which will eventually be taken to 1. We have also enforced that k is ultrasoft by neglecting it 
in front of M in the collinear gluon propagator. We shall evaluate Hfi.24(l in light cone coordinates. 
If we carry out first the integration over only the pole fc_ = kj_/k+ contributes. Then the only 
remaining singularities in the integrand are in the collinear gluon propagator. Hence, the absorptive 
piece can only come from its pole M 2 (l — z) — Mk + = kj_. If k + < M(l — z), then kj_ <~ M 2 (l — z) 
which implies fc_ ~ M. This contradicts the assumption that k is ultrasoft. Hence, \l\ must be 
expanded in the collinear gluon propagator. We then have 



1 




dk+8(M(\ -z)- k+) x 



K - E n + -f- + xl ^ . 



where we have introduced the change of variables |kjj = y/2k+x. Restricting ourselves to the 
ground state (n = 1) and using the techniques of reference |88| we obtain 

Im 



j o 

f°° / h h 

x ^ dx(2i; w (0)Is('f+x)-I 2 s ('f + : 




where 



7= } — z'=- =Ei- — -x X = E x =-— 6.27 

1 2 7 m 2 2N c C f m K ' 

This result can be recast in the factorized form given in |64| (eauation l6.2l) . 

Im (T$ So) (z)) = -vT' J dl + S s (l+)ImJ M (l + - M(l - z)) (6.28) 

ImJ M a+ - M(l - z)) = Tp (N 2 C - 1) ^S(M(1 - z) - 1+) (6.29) 
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We have thus obtained the 5- wave color octet shape function Ss{l+)- Analogously, for the P-wave 
color octet shape functions, we obtain from the second diagram of fig. 16.21 

Im (T0ff(z)y-r ] T'S jj 'Jdl + (s i fS Pi (l + ) + UV'- iaf) S P2 (l+)^ lmJ M (l+ - M(l - z)) 

(6-31) 

S P i{h) := £dx (Vo(0)/ P (-± + x)- IU 1 -J + *)) (6-32) 

Sp2{l+) ■= So dx {h+l x f (^ o(0) " 2 ^ io(0)/p( t +x)+ l2 ^i + x) ) (6 - 33) 

where 

= (2 " A) 4(i + »-) 3 ( 2(1 + = ' )(2 + 3 ' ! + (5 + + A) + 21-1 + A)2+ 

+ <rhy (<••' + -W - A2 » (- 1 + ( 1+ A i-)(/-A) + 

(6.34) 

Note that two shape functions are necessary for the P-wave case. 

The shape functions ((6.30(1 . 1(6.32(1 and (|6.33|l are ultraviolet (UV) divergent and require regu- 
larization and renormalization. In order to regulate them at this order it is enough to calculate the 
ultrasoft (US) loop (the integral over k in (|6.24|i 1 in P-dimensions, leaving the bound state dynamics 
in 3 space dimensions (P = 4 — 2e). In fact, the expressions 116.30(1 and H6.32[) implicitly assume that 
dimensional regularization (DR) is used, otherwise linearly divergent terms proportional to ip1 {0) 
would appear (which make (|6.3()|) and ((6.32(1 formally positive definite quantities). This procedure, 
to use DR for the US loop only, was the one initially employed in |42|. There the following steps 
were performed: 

• In order to isolate the 1/e poles, Is and Ip were expanded up to 0(1/ z' 2 ) (the expansion of 
these functions up to C(l/z' 4 ) can be found in equations 1(6.51(1 and 1(6.52(1 ). 

• The result was subtracted and added to the integrand of 1(6.30(1 - 1(6.32(1 (for ((6.33(1 this is not 
necessary since the only divergent piece is independent of Ip). The subtracted part makes the 
shape functions finite. The added part contains linear and logarithmic UV divergencies. 

• The remaining divergent integrals were dimensionally regularized by making the substitution 
dx — > dx(x//ji)~ £ . That produced the 1/e poles displayed in formulas (16) of ref. which 
were eventually subtracted (linear divergences are set to zero as usual in DR). 

That last point was motivated by the fact that x ~ (k_L being the transverse momentum of the 
US gluon) but differs from a standard MS scheme. 

As was already mentioned in |42j . this regularization and renormalization scheme is not the 
standard one in pNRQCD calculations. Later, in 03], a regularization-renormalization procedure 
closer to the standard one in pNRQCD was used; which is the one we will use here. That latter 
procedure consists in regularizing both the US loop and the potential loops (entering in the bound 
state dynamics) in DR; then US divergences are identified by taking the limit P — > 4 in the US 
loops while leaving the potential loops in P dimensions |89j ; potential divergencies are identified by 
taken P — > 4 in the potential loops once the US divergencies have been subtracted. It turns out 
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that all divergencies in Sp2 are US and all divergencies in Ss are potential. Spi contains both US 
and potential divergencies. The potential divergences related with the bound state dynamics can be 
isolated using the methods of ref. (HOj • Following this procedure we obtain the following expressions 
for the singular pieces 

c 2 7 a s (a u )^C 2 fa 2 (n v ) 
fc+)Uo - V L f — (1 - A ) (- 2 + A (2 In 2 + 1)) • 




M 2 



pC 





(6.35) 



9tt 2 N c U 6 ; v 6 

+ 2o^/_l_ ln f/|' )+ ... ) 



\^ 2 "r" m 



cn s (uuh 3 mC 2 fa 2 (ii v ) ( 31 19 1 

*M*+)Uo * h# + A(41n2 + 4) - A 2 (21n2 + ~ 



9TT 2 N r m 



,2 



3^ 2 7V C ^ 

For simplicity, we have set D = 4 everywhere except in the momentum integrals. fi p , according to 
H6.21|l , is given by (i p — y/m(M(l — z)/2 — E%). ^t c and /i pc are the subtraction points of the US and 
potential divergencies respectively. For comparison, let us mention that when we set /i c = M yT — z 
and [i pc = y/mnc, as we will do, we obtain exactly the same result as in the procedure used in 
ref. 02] for what the potential divergences is concerned 2 ; for the US divergences there is a factor 
In (2^7 J °f difference with respect to that former scheme. 

The renormalization of that expressions is not straightforward. We will assume that suitable 
operators exists which may absorb the 1/e poles so that an M S- like scheme makes sense to define the 
above expressions, and discuss in the following the origin of such operators. In order to understand 
the scale dependence of (|6.35|) - (|6.37|) it is important to notice that it appears because the term k 2 ^ 
in the collinear gluon propagator is neglected in l|6.24|l . It should then cancel with an IR divergence 
induced by keeping the term k 2 ^, which implies assuming a size M 2 (l — z) for it, and expanding 
the ultrasoft scales accordingly. We have checked that it does. However, this contribution cannot 
be computed reliably within pNRQCD (neither within NRQCD) because it implies that the k- 
component of the ultrasoft gluon is of order M, and hence it becomes collinear. A reliable calculation 
involves (at least) two steps within the EFT strategy. The first one is the matching calculation of the 
singlet electromagnetic current at higher orders both in a s and in (kj_/M) 2 and k+/M. The second 
is a one loop calculation with collinear gluons involving the higher order singlet currents. Figure 1^31 
shows the relevant diagrams which contribute to the IR behavior we are eventually looking for. We 
need NNLO in a s , but only LO in the (kj_/M) 2 and k+/M expansion. These diagrams are IR finite, 
but they induce, in the second step, the IR behavior which matches the UV of H6.35H - H6.37ll . The 
second step amounts to calculating the loops with collinear gluons and expanding smaller scales in 
the integrand. We have displayed in fig. 16.41 the two diagrams which provide the aforementioned IR 
divergences. For the UV divergences that do not depend on the bound state dynamics, we need the 
matching at LO in a s (last diagram in Fig. I6.3J1 but NLO in k + /M and (k_L/M) 2 . 

The above means that the scale dependence of the leading order contributions of the color- 
octet currents is of the same order as the NNLO contributions in a s of the color-singlet current, a 

2 Wc assume that the correlation of scales advocated in |1] (see 1241 for the implementation in our framework) must 
also be taken into account here. 
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Figure 6.4: Diagrams which induce an IR scale dependence which cancels against the UV one of the 
octet shape functions. 
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calculation which is not available. One might, alternatively, attempt to resum logs and use the NLO 
calculation |91| as the boundary condition. This log resummation is non-trivial. One must take into 
account the correlation of scales inherent to the non-relativistic system 0] , which in the framework 
of pNRQCD has been implemented in ^JESIE], and combine it with the resummation of Sudakov 
logs in the framework of SCET |t>lllt>2ll6*3"llf>4| (see also [EE])- Correlations within the various scales 
of SCET may start playing a role here as well In any case, it should be clear that by only 

resumming Sudakov logs, as it has been done so far |62| . one does not resum all the logs arising in 
the color octet contributions of heavy quarkonium, at least in the weak coupling regime. 

Keeping this in mind, we can proceed and write the renormalized expressions for the shape 
functions. These renormalized expressions, in an MS scheme, read 




Sf S ik + ) = ^^(^)'{ 2*io(0)( mjl^ ) ( / ( +X) - I - + 2Aln 2) 



7 k+ , 7 2 \ ( ha s N c \ ( [°°f ~% k. 
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In ref. |43| an additional subtraction related to linear divergencies was made. This subtraction was 
necessary in order to merge smoothly with the results in the central region. We will also need this 
subtraction when discussing the merging at LO in the following sections. We use 



00 1 7 
dx- — ► -2 
z' 



2 m V 2 



which differ from the MS scheme by the subtraction of the second term in the square brackets. In 
that other scheme (sub) the expressions for the shape functions read 



?subn„ \ qMS/i. \ 1 as 



Sf," 2 6 (fe+) = ^(fc+) (6.44) 

The validity of the formulas for the shape functions is limited by the perturbative treatment of 
the US gluons. The typical momentum of these gluons in light cone coordinates turns out to be: 



(fc+,fe ± ,fe_) = \M(l-z),^2M(l-z) ( M(1 2 Z) -E^j,M(l-z)-2EA (6.45) 

Note that the typical k± is not fixed by the bound state dynamics only but by a combination of the 
latter and the end-point kinematics. Hence, the calculation is reliable provided that k± > lGeV, 
which for the T(15) system means z < 0.92. 

6.2.4 Comparison with experiment 

We apply here the results in this section to the T(IS') system. There is good evidence that the 
T(15) state can be understood as a weak coupling (positronium like) bound state [BTl IBTfl 1701 ITTT 
E2 EH E3 EH EH IS! • Hence, ignoring O(Aqcd) in the shape functions, as we have done, should 
be a reasonable approximation. 

We will denote the contribution in the upper end-point region by . It is given by 

rlT e dT e dV e 

^ = ^ps + (6 46) 

az az dz 

where CS and CO stand for color singlet and color octet contributions respectively. The color 
singlet contribution is the expression with the Sudakov resummed coefficient (I6.6|) . The color octet 
contribution is given by 



dT co , s , s fWMa\ f 



a s (fi u ) a s (fih) — — t / C(x - z)S s +p(x)dx (6.47) 



where fi u is the US scale, that arises from the couplings of the US gluons (see below for the expression 
we use). C(x — z) contains the Sudakov resummations explained in l6.2.2P . The (tree level) matching 
coefficients (up to a global factor) and the various shape functions are encoded in S$+p(x), 



3 These matching coefficients, provided in reference 1621 . become imaginary for extremely small values of z — 1, a 
region where our results do not hold anyway. We have just cut-off this region in the convolutions. 
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0.75 0.8 0.85 0.9 0.95 

Figure 6.5: End-point contribution of the spectrum, dT e / dz, with the shape functions renormalized 
in an MS scheme (blue dashed line) and in the sub scheme (red solid line). The points are the 
CLEO data |S] 



(35pi(M(1 - z)) + S P2 (M(1 - z))) (6.48) 



The shape functions S s , S P1 and S P2 may become S^ s , S P [ S and Sgf or S s s ub , Spf and S a pf 
depending on the subtraction scheme employed. We will use the following values of the masses for 
the plots: m = 4.81 GeV and M = 9.46 GeV. The hard scale fih is set to fih = M. The soft 
scale fi s — mCfOi s is to be used for the a s participating in the bound state dynamics, we have 

a s (n s ) = 0.28. The ultrasoft scale [i u is set to [i u — ^2M{\ — z) (4f(l — z) — Ei), as discussed 
in the previous subsection. We have used the Mathematica package RunDec 94 to obtain the (one 
loop) values of a s at the different scales. 

In fjgure ET51 we plot the end-point contribution (|6.46|) with the shape functions renormalized in an 
MS scheme (blue dashed line) and in the sub scheme (red solid line) , together with the experimental 
data | 95 | (we have convoluted the theoretical curves with the experimental efficiency, the overall 
normalization of each curve is taken as a free parameter). We see that a very good description of 
data is achieved and that both schemes are equally good for the description of the shape of the 
experimental data in the end-point region. This nice agreement with data is an encouraging result. 
But still remains to be seen if it is possible to combine these results, for the end-point region, with 
the ones for the central region (where the NRQCD formalism is expected to work) . This will be the 
subject of section 1531 




70 



Chapter 6. Radiative heavy quarkonium decays 




Figure 6.6: Diagrams which require a Vi(^S\) operator for renormalization. The solid circle stands 
for either the Os^Sq) or Os( 3 Pj) operator, the crossed box for either the chromomagnetic (<■ ) or 
chromoelectric (O ) interaction in fig. 16.21 the empty box for the octet Coulomb potential, and the 
thin solid lines for free QQ propagators. 



6.2.5 Calculation of T(1S) NRQCD color octet matrix elements 

The calculation of the shape functions can be easily taken over to provide a calculation of (T(15)| 
8 ( 1 So)|T(lS')) and (T(1S)\0 8 ( 3 Pj)\T(1S)) assuming that ma 2 > A QCD is a reasonable approx- 
imation for this system. Indeed, we only have to drop the delta function (which requires a further 
integration over fc + ) and arrange for the suitable factors in (|6.28(1 and H6.31[l . We obtain 

/•OO 

<T(lS)|0 8 ( 1 ^o)|T(l^)) - -2T 2 {N 2 C - 1) / dk + S s {k+) 

Jo 

(T(1S)\0 8 ( 3 Pj)\T(1S)) = _ < 2J + 1 ) T f(^ - 1) J°° dk+Spi(k+) ((U9) 
where we have used 

poo 2 r°° 

/ dk + S P2 (k + ) = - dk + S P1 (k+) (6.50) 
Jo Jo 

The expressions above contain UV divergences which may be regulated by calculating the ultrasoft 
loop in D dimensions. These divergences can be traced back to the diagrams in fig. 16. 61 and fig. 16.71 
Indeed, if we expand Is and Ip for z' large, we obtain 



AV 



+ \ -\ + A(21n2 + 1) + A 2 (-41n2) + -C(3)A J J + G(-^)| (6.51) 

'-i/?So-«{^(-!^HH)^ 

1 - A+ — (-6 + tt 2 )A 2 ) -L + - (-5 + A + A 2 (2-81n2) + 8Aln2+ 



+A 3 (-41n2 + 3C(3)))i T + ^)} 



(6.52) 



It is easy to see that only powers of 1/z' up to order four may give rise to divergences. Moreover, 
each power of 1/z' corresponds to one Coulomb exchange. Taking into account the result of the 
following integral, 



i _ iW7 2 \, 2 ~ 2£ r 2 (i- £ ) 



dk + I dx(2k + x)-"— = 2 1 -^[^-) ^ ^ - + 2e-2j (6.53) 



o 
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Figure 6.7: Diagrams which require a Oi( 3 Si) operator for renormalization. Symbols are as in fig. 

una 



we see that only the 1/ z' 2 and 1/z' 4 terms produce divergences. The former correspond to diagrams 



in fig. 16.61 and the latter to fig. 16.71 which can be renormalized by the operators "Pi( 3 5i) and 
0i( 3 Si) respectively. It is again important to notice that these divergences are a combined effect 
of the ultrasoft loop and quantum mechanics perturbation theory [potential loops |96p and hence 
it may not be clear at first sight if they must be understood as ultrasoft (producing log/x„ in the 
notation of refs. |1U1 1921 HQ ) or potential (producing \ogfi p in the notation of refs. |101 1921 HQ ). 
In any case, the logarithms they produce depend on the regularization and renormalization scheme 
used for both ultrasoft and potential loops. Remember that the scheme we use here is not the 
standard one in pNRQCD [HH1 E3 HI3 E2 CHI EHj • In the standard scheme the ultrasoft divergences 
(anomalous dimensions) are identified by dimensionally regulating both ultrasoft and potential loops 
and subsequently taking D — > 4 in the ultrasoft loop divergences only. If we did this in the present 
calculation we would obtain no ultrasoft divergence. Hence, in the standard scheme there would be 
contributions to the potential anomalous dimensions only. The singular pieces in our scheme are 
displayed below 
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(6.54) 



12 9 

With these expressions we obtain the following estimates for the value of the matrix elements 

<T(15)|0 8 ( 1 5o)|T(I5))| ai=m ~ 0.004Gey 3 (6.55) 
(T(15)|O 8 ( 3 F )|T(15))| Ai=M ~ 0.08 GeV 5 (6.56) 
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remember that the above numbers are obtained in an MS scheme from dimensionally regularized 
US loops only. The value we assign to the iS-wave matrix element is compatible with the recent 
(quenched) lattice determination (hybrid algorithm) |9"%] . 

6.3 Merging the various contributions to the spectrum 

Now, with the, for a long time elusive, end-point region of the spectrum well described, it is the 
time to put together all the contributions to the spectrum and see if a good description of data is 
achieved. As was already explained, the contributions to the spectrum can be split into direct ( dlr ) 
and fragmentation (f ra 9) 



The fragmentation contributions arc those in which the photon is emitted from the decay products of 
the heavy quark (final state light quarks), these contributions where first taken into account in |82| 
and further studied in 57 ; while the direct contributions are those in which the photon is emitted 
directly from the heavy quark itself. Although this direct-fragmentation splitting is correct at the 
order we are working it should be refined at higher orders. We discuss each of these contributions 
in turn, in the two following subsections. 

6.3.1 Direct contributions 

The approximations required to calculate the QCD formula (|6.15|) are different in the lower end-point 
region (z — > 0), in the central region (z ~ 0.5) and in the upper end-point region (z — > 1). 

In the lower end-point region the emitted low energy photon can only produce transitions within 
the non-relativistic bound state without destroying it. Hence the direct low energy photon emission 
takes place in two steps: (i) the photon is emitted (dominantly by dipole electric and magnetic 
transitions) and (ii) the remaining (off-shell) bound state is annihilated into light hadrons. This 
lower end-point contribution goes to zero, for z — > 0, as z 3 , while the leading order NRQCD result 
goes to zero as z (see HH1 llUUj for a recent analysis of this lower end-point region in QED). As 
was already mentioned, at some point the direct photon emission is overtaken by the fragmentation 
contributions [821 157|. In practice this happens about z ~ 0.4, namely much before than the z 3 
behavior of the low energy direct photon emission can be observed, and hence we shall neglect the 
latter in the following. 

For z away from the lower and upper end-points (0 and 1 respectively), no further scale is 
introduced beyond those inherent of the non-relativistic system. The integration of the scale m 
in the time ordered product of currents in 1)6. 15[) leads to local NRQCD operators with matching 
coefficients which depend on m and z. We will summarize here the known results for the central 
region (we denote the direct contributions in the central region by T c ). At leading order one obtains 



dT dT dir dTf ra 3 



(6.57) 



dz dz dz 
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2-z z(l-z) 
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ln(l - z), 



(6.58) 



where 4 




(6.59) 



4 Note added: The (^r) factor was missing in formula (4) of 1441 (and in a previous version of the thesis). This 
was just a typo, not affecting any of the subsequent results. We thank A.Vairo for help in identifying it. 
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The a s correction to this rate was calculated numerically in ref. |91| . The expression corresponding 
to (|6.59() in pNRQCD is obtained at lowest order in any of the possible regimes by just making the 
substitution 

(V Q (nS)\0 1 ( 3 S 1 )\V Q (nS)) = 2N c \^ n0 (0)\ 2 , (6.60) 

where ip n o(0) is the wave function at the origin. The final result coincides with the one of the early 
QCD calculations [461 147| . We will take the Coulomb form V'io(O) = j 3 /tt for the LO analysis of 
T(15). 

The NLO contribution in the original NRQCD counting [3] is v 2 suppressed with respect to 
It reads 

dz ivy TO 4 

In the original NRQCD counting or in the weak coupling regime of pNRQCD the new matrix 
element above can be written in terms of the original one jl()l| 5 

(V Q (nS)\V i eS 1 )\V Q (nS)) _ / M - 2m \ {V Q {nS)\Q^ S 1 )\V Q (nS)) , | Q , 2)) 
m 4 \ m J m 2 

The matching coefficient can be extracted from an early calculation j^S] (see also |54| ^ . It reads 

16 2 2 (- , , 1 
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(6.64) 



The contributions of color octet operators start at order v A . Furthermore, away of the upper 
end-point region, the lowest order color octet contribution identically vanishes [Uj. Hence there is 
no l/a s enhancement in the central region and we can safely neglect these contributions here. 

If we use the counting a s ((J,h) ~ v 2 , a s (/i s ) ~ v (remember that /j,h ~ m and jj, s ~ mv are the 
hard and the soft scales respectively) for the T(15*), the complete result up to NLO (including v 2 
suppressed contributions) can be written as 

^El = ^lo + d^NLO + dT lo, aB ( g g 5) 

dz dz dz dz 

The first term consist of the expression (16.58(1 with the Coulomb wave function at the origin l|6.60|) 
including corrections up to O (a s ([i s )) 2 [11)21 11U3] . the second term is given in (|6.61l) . and the 
third term consists of the radiative O (a s (fih)) corrections to l|6.58[l which have been calculated 

5 In the strong coupling regime of pNRQCD an additional contribution appears 1871 
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numerically in |91j . Let us mention at this point that the O (a s (/x s )) 2 corrections to the wave 
function at the origin turn out to be as large as the leading order term. This will be important for 
the final comparison with data at the end of the section. Note that the standard NRQCD counting 
we use does not coincide with the usual counting of pNRQCD in weak coupling calculations, where 
ots{Hh) ~ &s([i-s) ~ a s (mv 2 ). The latter is necessary in order to get factorization scale independent 
results beyond NNLO for the spectrum and beyond NLO for creation and annihilation currents. 
However, for the T(15) system (and the remaining heavy quarkonium states) the ultrasoft scale mv 2 
is rather low, which suggests that perturbation theory should better be avoided at this scale |72|. 
This leads us to standard NRQCD counting. The factorization scale dependences that this counting 
induces can in principle be avoided using renormalization group techniques 0] HOI 1241 1921 125] . In 
practice, however, only partial NNLL results exists for the creation and annihilation currents |26II27| 
(see |1U4| for the complete NLL results) , which would fix the scale dependence of the wave function 
at the origin at 0(a 2 (mv)). We will not use them and will just set the factorization scale to m. 

The upper end-point region of the spectrum has been discussed in great detail in the previous 
section. As we have seen there, different approximations, with respect to the ones for the central 
region, are needed here. It is not, by any way, obvious how the results for the central and for 
the upper end-point regions must be combined in order to get a reliable description of the whole 
spectrum. When the results of the central region are used in the upper end-point region, one misses 
certain Sudakov and Coulomb resummations which are necessary because the softer scales Myl — z 
and M(l — z) become relevant. Conversely, when results for the upper end-point region are used in 
the central region, one misses non-trivial functions of z, which are approximated by their end-point 
[z ~ 1) behavior. We will explain, in the remaining of this subsection, how to merge these two 
contributions. 



Merging the central and upper end-point regions 

One way to proceed with the merging is the following. If we assume that the expressions for the 
end-point contain the ones of the central region up to a certain order in (1 — z), we could just 
subtract from the expressions in the central region the behavior when z — > 1 at the desired order 
and add the expressions in the end-point region. Indeed, when z — > 1 this procedure would improve 
on the central region expressions up to a given order in (1 — z), and when z belongs to the central 
region, they would reduce to the central region expressions up to higher orders in a s . This method 
was used in ref. |64| and in ref. |43) . In ref. |64) only color singlet contributions were considered 
and the end-point expressions trivially contained the central region expressions in the limit z — > 1. 
In ref. |43| color octet contributions were included, which contain terms proportional to (1 — z). 
Hence, the following formula was used 

1 dT dlr 1 dT c rn ( 1 dT e nq \ ( 1 dY e rn , , , NW ,\ , 

^ ' 2£ - z + -fQ- - z (4 + 2 log (1 - z)) (1 - z)\ (6.66) 



Tq dz To dz \Tq dz J \To dz 

Even though a remarkable description of data was achieved with this formula (upon using a suitable 
subtraction scheme, the sub scheme described in the previous subsection), this method suffers from 
the following shortcoming. The hypothesis that the expressions for the end-point contain the ones 
for the central region up to a given order in (1 — z) is in general not fulfilled. As we will see below, 
typically, they only contain part of the expressions for the central region. This is due to the fact 
that some a s (nh) in the central region may soften as a s (M(l — z)), others as a s (M y/l — z) and 
others may stay at a s (nh) when approaching the end-point region. In a LO approximation at the 
end-point region, only the terms with the a s at low scales would be kept and the rest neglected, 
producing the above mentioned mismatch. We shall not pursue this procedure any further. 
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Let us look for an alternative. Recall first that the expressions we have obtained for the upper 
end-point region are non-trivial functions of M(l — z), Myl — z, ma s (mv) and ma^mu), which 
involve a s at all these scales. They take into account both Sudakov and Coulomb resummations. 
When z approaches the central region, we can expand them in a s (M\/l — z), a s (M(l — z)) and the 
ratio ma s (mv)/M^/l — z. They should reduce to the form of the expressions for the central region, 
since we are just undoing the Sudakov and (part of) the Coulomb resummations. Indeed, we obtain 
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f(x,z) = (1 - a;) log 



(1 - z) log 



(6.68) 



(6.69) 
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M 2 (l-z) 2 , 
-64C/(2 — A)/9 in the sub scheme). 



M 2 {l-x) 2 ^ 

-N c - 136C/(2- A)/9 (in an MS scheme; it becomes A ■ 
In the next paragraph we explain how to obtain these expressions for the shape functions in the 
central region. 

First consider the S'-wave octet shape function 



hfe+x):= I d a x^ 10 (x)( 1- 



(6.70) 



x,0 



p 2 /m + V , V 



a s /(2AT c |r|). When z approaches the central region, k + ~ M(l — z) ^> —E-y 

the typical three momentum in the bound 



and the larger three momentum scale is Mi/1 — z 3> 7, 
state. Therefore we can treat the Coulomb potential in H6.70|l as a perturbation when it is dominated 
by this scale. It is convenient to proceed in two steps. First we write h a = h s + (V a — V s ), where 
h s = p 2 /m + V s , V s = — a s C//|r|, and expand V — V s . This allows to set h s — E\ to zero in 
the left-most propagator and makes explicit the cancellation between the first term in the series 
and the first term in (|6.70() . It also makes explicit that the leading term will be proportional to 
a s (M v / l — z). Second, we expand V s in h s = p 2 /m + V s . In addition, since My/1 — z ^> 7, the 
wave function can be expanded about the origin. Only the first term in both expansion is relevant 
in order to get H6.68JI . Consider next the P-wave shape functions 



Ipfe+x) '=-1 I d a xxV 10 (x) 



1 - 



h -E 1 + ^ 



(6.71) 
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In order to proceed analogously to the S- wave case, we have first to move the x ! away from the wave 
function ^ 

M% +x)= Vio(O) + [ d 3 x^ 10 (x) ' 



Ei 



xV 



E 1 + 



2V J 



k+ 

2 



(6.72) 



For the left-most propagators we can now proceed as before, namely expanding V — V s . Note that 
the leading contribution in this expansion of the second term above exactly cancels against the first 
term. Of the remaining contributions of the second term only the next-to-leading one (0(a s )) is 
relevant to obtain (|6.68|) . Consider next the leading order contribution in this expansion of the last 
term. It reads 



3 m 



d 3 xV>io(x) I V 
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(6.73) 



Now we proceed as before with the left-most propagators, namely expanding V — V s . The leading 
order contribution of the first term above produces the relativistic correction 0(v 2 ) of (|6.68|) . The 
next-to-leading contribution of this term and the leading order one of the second term are 0(a a ) 
and also relevant to (|6.68() . The next-to- leading order contribution of the last term in (|6.72() in the 
V ~ V s expansion of the left-most propagator is also 0(a s ) and relevant to Ij6.68|l . 

Returning now to equations l|6.67(l - H6.68(l . we see that the color singlet contribution reproduces 
the full LO expression for the central region in the limit z — > 1. The color octet shape functions 
Spi and Sp2 give contributions to the relativistic corrections l|6.61|l . and Sp2 to terms proportional 
to (1 — z) in the limit z — > 1 of (|6.58l) as well. We have checked that, in the z — » 1 limit, both 
the (1 — z) ln(l — z) of l|6.58|l and the ln(l — z) of the relativistic correction (|6.61() are correctly 
reproduced if \jl c ~ My/1 — z, as it should. All the color octet shape functions contribute to the 
0{a s {fJbh)) correction in the first line of (|6.68(l . There are additional 0{a s {nh)) contributions coming 
from the expansion of the (Sudakov) resummed matching coefficients of the color singlet contribution 
and of the Sp2 color octet shape function. The a s log(l — z) in l|6.67|l reproduces the logarithm in 
dT c LO , a Jdz. 

We propose the following formula 
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(6.74) 



This formula reduces to the NRQCD expression in the central region. When we approach the upper 
end-point region the second terms in each of the parentheses are expected to cancel corresponding 
terms in the z-* 1 limit of the expression for the central region up to higher order terms (in the 
end-point region counting). Thus, we are left with the resummed expressions for the end-point (up 
to higher order terms). 

There are of course other possibilities for the merging. For instance, one may choose a z\ below 
which one trusts the calculation for the central region and a z 2 above which one trusts the end-point 
region calculation, and use some sort of interpolation between z\ and z 2 (see for instance |105p . This 
would have the advantage of keeping the right approximation below z\ and beyond z 2 unpolluted, at 
the expense of introducing further theoretical ambiguities due to the choice of z\ and z 2l and, more 
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important, due to the choice of the interpolation between Z\ and z%. We believe that our formula 
(|6.74|l is superior because it does not introduce the above mentioned theoretical ambiguities. The 
price to be paid is that the expressions from the central region have an influence in the end-point 
region and vice- versa. This influence can always be chosen to be parametrically subleading but large 
numerical factors may make it noticeable in some cases, as we shall see below. 

Merging at LO 

If we wish to use only the LO expressions for the central region, we should take (|B~H7|l and ijfTBgj) at 

LO, namely 

r z ( 2 - Uog (^))) {1 - z) (6J5) 

and substitute them in (|6.74|) . Unexpectedly, the results obtained with this formula in the central 
region deviate considerably from those obtained with formula (|6.58() (see fig. 16.8(1 . This can be 
traced back to the fact that the a s y/l — z corrections in l(6.68|l are enhanced by large numerical 
factors, which indicates that the merging should better be done including a s (p,h) corrections in the 
central region, as we discuss next. Alternatively, we may change our subtraction scheme in order 
to (partially) get rid of these contributions. With the new subtraction scheme (sub), described in 
the preceding section, the situation improves, although it does not become fully satisfactory (see fig. 
16. 8f) . This is due to the fact that some a s ^/l — z terms remain, which do not seem to be associated 
to the freedom of choosing a particular subtraction scheme. In spite of this the description of data 
turns out to be extremely good. In figure l6~$l we plot, using the sub scheme, the merging at LO 
(solid red line) and also, for comparison, equation (|6.66|l (blue dashed line). We have convoluted 
the theoretical curves with the experimental efficiency and the overall normalization is taken as a 
free parameter. 

Merging at NLO 

If we wish to use the NLO expressions for the central region ((6.65(1 . we should take all the terms 
displayed in l|6.67|l - l|6.68|l and substitute them in (|6.74|) . Unlike in the LO case, for values of z in 
the central region the curve obtained from (|6.74(l now approaches smoothly the expressions for the 
central region 1(6.65(1 as it should. This is so no matter if we include the a^(^i s ) corrections to the 
wave function at the origin in dT c LO /dz, as we in principle should, or not (see figs. 16.101 and I6.11JI . 
However, since the above corrections are very large, the behavior of the curve for z — > 1, strongly 
depends on whether we include them or not (see again figs. 16.101 and I6.11|) . We believe that the 
two possibilities are legitimate. If one interpretes the large a|?(/i s ) corrections as a sign that the 
asymptotic series starts exploding, one should better stay at LO (or including a s (p s ) corrections). 
However, if one believes that the large a^(iJ, s ) corrections are an accident and that the a^(/j, s ) ones 
(see EO] for partial results) will again be small, one should use these a^(fi s ) corrections. We 
consider below the two cases. 

If we stay at LO (or including a s (p s ) corrections) for the wave function at the origin, the 
curve we obtain for z — > 1 differs considerably from the expressions for the end-point region 1(6.4611 
(see fig. 16.10(1 . This can be traced back to the a s ^/l — z term in 1(6.68(1 again. This term is 
parametrically suppressed in the end-point region, but, since it is multiplied by a large numerical 
factor, its contribution turns out to be overwhelming. This term might (largely) cancel out against 
higher order contributions in the end-point region, in particular against certain parts of the NLO 
expressions for the color singlet contributions, which are unknown at the moment. 

If we use the wave function at the origin with the al(^s) corrections included, the curves we 
obtain for z — > 1 become much closer to the expressions for the end-point region 1(6.46(1 (sec fig 
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Figure 6.8: Merging at LO. The solid red line is the NRQCD expression (|6.58[) . The dot-dashed 
curves are obtained using an MS scheme: the pink (light) curve is the end-point contribution l|6.46|) 
and the black (dark) curve is the LO merging. The dashed curves are obtained using the sub 
scheme: the green (light) curve is the end-point contribution H6.46J) and the blue (dark) curve is the 
LO merging. 




Figure 6.9: Direct contribution to the spectrum. The solid red line corresponds to the LO merging 
and the blue dashed line corresponds to equation H6.66fl . The points are the CLEO data |95j . 
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Figure 6.10: Merging at NLO (using an MS scheme and the wave function at the origin at LO). The 
solid red line is the NRQCD result (|6.65|1 . the blue (light) dashed curve is the end-point contribution 
(|6.46J) and the black (dark) dashed curve is the NLO merging. 

I6.11[l . Hence, a good description of data is obtained with no need of additional subtractions 6 , as 
shown in figure 15". 121 (as usual experimental efficiency has been taken into account and the overall 
normalization is a free parameter). This are now good news. Because this final curve incorporates 
all the terms that are supposed to be there according to the power counting. 

6.3.2 Fragmentation contributions 

The fragmentation contributions can be written as 

dTf ra 3 K f 1 dx 



£ / ^C a (x)D a ,(^M), (6.76) 



dz 

a=q,q,g 

where C a represents the partonic kernels and D aj represents the fragmentation functions. The 
partonic kernels can again be expanded in powers of v |57| 

C a =Y J C a [Q] (6.77) 
Q 

The leading order term in v is the color singlet rate to produce three gluons 



c, [ftfaa = g«? + + 2 l -± ln(1 _ „ _ 2 jizii_ ln(1 _ „ 



2 



e One might be worried about the big difference that the corrections to the wave function at the origin introduce in 
the result. In that sense let us mention that when we analyze the electromagnetic decay width of T(1S) (T(T(1S) — > 
e+ e~ ), the formulas needed to compute the width can be found, for instance, in HUtil ). with the same power counting 
we have employed here, the result we obtain is 5.24 ■ 10" 7 GeV if we do not include the Qj(/i s ) corrections and 
1.17 • 10 _6 GeV if we do include them. This is to be compared with the experimental result 1.32 • 10 — 6 GeV I1U7I . 
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Figure 6.11: Merging at NLO (using an MS scheme and the wave function at the origin with the 
a%(li s ) corrections included). The solid red line is the NRQCD result 1)6. 65[) . the blue (light) dashed 
curve is the end-point contribution 1)6. 46[) and the black (dark) dashed curve is the NLO merging. 




Figure 6.12: Direct contribution to the spectrum using the NLO merging (in an MS scheme and 
the wave function at the origin with the «g(/i s ) corrections included). The points are the CLEO 
data U2J. 
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m 2 

The color octet contributions start at order v A but have a enhancement with respect to <|6.78[1 

C g [O 8 CS 0) ] = ^6(1 - z) (y Q (nSWS a) \V Q (nS)) 
1 J 3 m 

C g [WPj)] = ^5(1 - z) (y Q (nS)\Os(y\V Q (nS)) 

3 777- 

C q WSl )} = ^5(1 - z) (y Q (nS)\Os(y\V Q (nS)) 
1 J 3 m z 

Then the color singlet fragmentation contribution is of order a 3 D g ^ and the color octet frag- 
mentation are of order v^ct^Dg^ ( 1 S'o and 3 Pj contributions) or v^a^Dg^ ( 3 Si contribution). 
We can use, as before, the counting v 2 ~ a s to compare the relative importance of the different 
contributions together with the existing models for the fragmentation functions |l()8j . The latter 
tell us that D q ^ 1 is much larger than D g ^ 7 . This causes the 0(v 4 a 2 D q ^) 3 Si octet contribution 
to dominate in front of the singlet 0{a 3 D g ^ 1 ) and the octet 0(v 4 a 2 D g ^ 7 ) contributions. In fact, 
a s D 9 _> 7 is still larger than D g ^ 11 so we will include in our plots the a s corrections to the color 
octet contributions (|fi.79|) proportional to -D q ^ 7 , which have been calculated in |S7]. In addition, 
the coefficients for the octet 3 Pj contributions have large numerical factors, causing these terms to 
be more important than the color singlet contributions. Let us finally notice that the a s corrections 
to the singlet rate will produce terms of 0{a s D q ^^), which from the considerations above are ex- 
pected to be as important as the octet 3 Si contribution. These a s corrections to the singlet rate 
are unknown, which results in a large theoretical uncertainty in the fragmentation contributions. 

For the quark fragmentation function we will use the LEP measurement 



D qj (z,fi) 
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(6.80) 



Ml, . D t ^_l + (l-z) 2 _ „ _ n14+0 .43 



C = -l-ln(-f) ; P qi (z) = i '- ; ^ = 0.14™GeV (6.81) 

and for the gluon fragmentation function the model |110j . These are the same choices as in |64) . 
However, for the O^CSq) and Os( 3 -Po) matrix elements we will use our estimates H6.55M6.5fi0 . 
Notice that we do not assume that a suitable combination of these matrix elements is small, as it 
was done in [Hlj. The 08( 3 <Si) matrix element can be extracted from a lattice determination of the 
reference |98| . Using the wave function at the origin with the a 2 (/j, s ) corrections included, we obtain 
(we use the numbers of the hybrid algorithm) , 

<T(1S*)|0 8 ( 3 5i)|T(15))| ai=m ~ 0.00026 GeV 3 (6.82) 

which differs from the estimate using NRQCD v scaling by more than two orders of magnitude: 

{T(lS)\O s ( 3 S 1 )\T(lS))\^ M ~ v 4 (T(15)|0 1 ( 3 5 1 )|T(15))| ai=m ~ 0.02G e y 3 (6.83) 

(we have taken v 2 ~ 0.08), which was used in ref. [Hi]- The description of data turns out to be 
better with the estimate (|6.83fl . However, this is not very significant, since, as mentioned before, 
unknown NLO contributions are expected to be sizable. 

In the z — > region soft radiation becomes dominant and the fragmentation contributions com- 
pletely dominate the spectrum in contrast with the direct contributions |82|. Note that, since the 



82 



Chapter 6. Radiative heavy quarkonium decays 



fragmentation contributions have an associated bremsstrahlung spectrum, they can not be safely 
integrated down to z — 0; that is f dz dr dz 3 is not an infrared safe observable. In any case we 
are not interested in regularizing such divergence because the resolution of the detector works as a 
physical cut-off. 




6.3.3 The complete photon spectrum 

We can now compare the theoretical expressions with data in the full range of z. First note that 
formula H6.74[l requires dT e j dz for all values of z. The color octet shape functions, however, were 
calculated in the end-point region under the assumption that Myl-z ~ 7, and the scale of the a s 
was set accordingly. When z approaches the central region Myl — z S> 7, and hence some a s will 
depend on the scale Ma/1 — z and others on 7 (we leave aside the global a s (fj, u )). In order to decide 
the scale we set for each a s let us have a closer look at the formula (|6.68|1 . We see that all terms 
have a common factor 7 s . This indicates that one should extract 7 3 factors in the shape functions, 
the a s of which should stay at the scale /i s . This is achieved by extracting 7 s / 2 in Is and Ip. If 
we set the remaining a s to the scale \i v = y/m{M{l — z)/2 — E\), we will reproduce l|6.68[) when 
approaching to the central region, except for the relativistic correction, the a s of which will be at 
the scale n P instead of at the right scale fx s . We correct for this by making the following substitution 

- (6.84) 

Notice that the replacements above are irrelevant as far as the end-point region is concerned, but 
important for the shape functions to actually (numerically) approach the expressions (|6.68l) in the 
central region, as they should. 

The comparison with the experiment is shown in figures 15 . 1 31 and 16 . 1 41 These plots are obtained 
by using the merging formula l|6.74|l at NLO with the a 2 (/i s ) corrections to the wave function at the 
origin included for the direct contributions plus the fragmentation contributions in subsection !6.3.2l 
including the first a s corrections in C q and using the estimate l|6.83[l for the (T(lS')|C'8( 3 S'i)|T(15)) 
matrix element. The error band is obtained by replacing /i c by v / 2 ±1 /x c . Errors associated to 
the large a 2 (/i s ) corrections to the wave function at the origin, to possible large NLO color singlet 
contributions in the end-point region and to the fragmentation contributions are difficult to estimate 
and not displayed (see the corresponding sections in the text for discussions). The remaining error 
sources are negligible. In figure 16". 131 as usual, experimental efficiency has been taken into account 
and the overall normalization is a free parameter. Figure 16.141 compares our results with the new 
(and very precise) data from CLEO 0H|- This plot takes into account the experimental efficiency 
and also the resolution of the experiment (the overall normalization is a free parameter). 

We can see from figures l6*.13l and l6.14l that. when we put together the available theoretical results, 
an excellent description of data is achieved for the whole part of the spectrum where experimental 
errors are reasonable small (recall that the error bars showed in the plots only take into account the 
statistical errors, and not the systematic ones |§S10Hj). Clearly then, our results indicate that the 
introduction of a finite gluon mass |1 1 1 j is unnecessary. One should keep in mind, however, that in 
order to have the theoretical errors under control higher order calculations are necessary both in the 
direct (end-point) and fragmentation contributions. 

Let us mention that the inclusion of color octet contributions in the end-point region, together 
with the merging with the central region expression explained here, may be useful for production 
processes like inclusive J/tp production in e+e~ machines [1121 11051 1113| . 
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Figure 6.13: Photon spectrum. The points are the CLEO data [HH]- The solid lines are the NLO 
merging plus the fragmentation contributions: the red (light) line and the blue (dark) line are 
obtained by using fgS)) and for (T(l^|0 8 ( 3 Si)|T(LS)) respectively. The grey shaded region 

is obtained by varying [i c by v / 2 ±1 /i c . The green shaded region on the right shows the zone where 
the calculation of the shape functions is not reliable (see subsection I6.2.3|) . The pink dashed line is 
the result in [HI], where only color singlet contributions were included in the direct contributions. 
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Figure 6.14: Photon spectrum. The points are the new CLEO data The red solid line is the 
NLO merging plus the fragmentation contributions, using (|6~%5|) for (T(lS , )|C 8 ( 3 S'i)|T(lS')). The 
grey shaded region is obtained by varying fx c by v / 2 ±1 /x c . The green shaded region on the right 
shows the zone where the calculation of the shape functions is not reliable (see subsection 16 . 2 . 3|) . 
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6.4 Identifying the nature of heavy quarkonium 

As we have just seen in the previous section, the photon spectrum in the radiative decay of the T(1S) 
can be well explained theoretically. This fact, together with the recent appearance of measurements 
of the photon spectra for the T(2S) and T(3S) states motivates us to try to use these radiative 
decays to uncover the properties of the decaying heavy quarkonia. 

As it has already been explained, the interplay of Aqcd with the scales mv and mv 2 dictates 
the degrees of freedom of pNRQCD. Two regimes have been identified: the weak coupling regime, 
Aqcd < mv 2 , and the strong coupling regime, mv 2 <C Aqcd ^ mv. Due to the fact that none of the 
scales involved in these hierarchies are directly accessible experimentally, given a heavy quarkonium 
state, it is not obvious to which regime it must be assigned. Only the T(15) appears to belong to 
the weak coupling regime, since weak coupling calculations in a s (mv) converge reasonably well. The 
fact that the spectrum of excitations is not Coulombic suggests that the higher excitations are not 
in the weak coupling regime, which can be understood from the fact that O(Aqcd) effects in this 
regime are proportional to a high power of the principal quantum number |114Mll5) . Nevertheless, 
there have been claims in the literature, using renormalon-based approaches, that also T(25*) and 
even T(35) can also be understood within the weak coupling regime |78l 1741 ITU)] . We will see that 
the photon spectra in semi-inclusive radiative decays of heavy quarkonia to light hadrons provide 
important information which may eventually settle this question. 

We start by writing the radiative decay rate for a state with generic principal quantum number 
n. Again we split the decay rate into direct and fragmentation contributions 

din = rfi^ + 

dz dz dz 

here z = 2E 1 /M n (M n is the mass of the heavy quarkonium state). We shall now restrict our 
discussion to z in the central region, in which no further scale is introduced beyond those inherent 
of the non-relativistic system. We write the spectrum in the following compact form 

' /ir =E C [2]W% (6-86) 



dz ^ m° a 

a=q,q : g Jz Q Q 



(6.87) 



where Q is a local NRQCD operator, 6q is an integer which follows from the dimension of Q and 
(Q)n '■— {VQ(nS)\Q\VQ(nS)}. It is important for what follows that the Jq(z) are universal and do 
not depend on the specific bound state n. Due to the behavior of the fragmentation functions above, 
the fragmentation contributions are expected to dominate the spectrum in the lower z region and 
to be negligible in the upper z one. In the central region, in which we will focus on, they can always 
be treated as a perturbation, as we will show below. 

Let us first consider the weak coupling regime, for which the original NRQCD velocity counting 
holds 3 (this is the situation described in the previous sections of this chapter, we recall here 
some of the arguments for an easier reading). The direct contributions are given at leading order 
by the 0\ ( 3 Si) operator; the next-to-leading order (NLO) (v 2 suppressed) term is given by the 
T 3 ! ( 3 Si) operator. The contributions of color octet operators start at order v 4 and are not a~ 1 (m) 
enhanced in the central region. The fragmentation contributions are more difficult to organize since 
the importance of each term is not only fixed by the velocity counting alone but also involves the 
size of the fragmentation functions. It will be enough for us to restrict ourselves to the LO operators 
both in the singlet and octet sectors. The LO color singlet operator is 0\ ( 3 5i) as well. The leading 
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color octet contributions are v 4 suppressed but do have a a^^m) 
to the singlet ones here. They involve Og ( 3 <Si)j Cs (^o) and Og 
the NRQCD expression (at the order described above) reads 



~ 1 jv 2 enhancement with respect 
3 Pq) . Then in the central region, 



C[ [ 3 S 1 ] (z) 



, (0 8 (^q))» , (08( 3 P O )) T 
-/o 8 (iSo)( z ) + /Os( 3 P.7)( 2 ) "4 



(6.88) 



If we are in the strong coupling regime and use the so called conservative counting, the color octet 
matrix elements are suppressed by v 2 rather than by w 4 . Hence we should include the color octet 
operators in the direct contributions as well. In practise, this only amounts to the addition of CVs 
to the /o 8 's. Furthermore, /o 1 (8 Sl )(z), fo s (^S )(z) and fo B (spj)(z) are proportional to D gl (x,m), 
which is small (in the central region) according to the widely accepted model |110| . /o 8 ( 3 s 1 )( 2; ) is 
proportional to D ql (x,m), which has been measured at LEP |109| . It turns out that numerically 
/o 8 ( 3 Si) ( z ) ~ Cs[ 3 Si](z) in the central region. Therefore, all the LO fragmentation contributions 
can be treated as a perturbation. Consequently, the ratio of decay widths of two states with different 
principal quantum numbers is given at NLO by 



dz 
dT r 

dz 



(Oi( 3 5i))„ L , C[ [ 3 S X ] (z)^ (3Sl) 



(Oi( 3 5i)) r I C 1 [ 3 S 1 ](z) m 2 



~' < -0 8 ( 3 S 1 )+ ^ - l< - 



Ci [ 3 Si] (z) 



Ci [ 3 S 1 ] (z) 



CM 1 So)" 



where 



fo a PPj){z) n 



O 8 ( a P ) 



d [ 3 5i] (z) m 2 



(Q)n 



(Qh 



(0 1 (3S 1 )) n (0 1 ( 3 S 1 )) r 



(6.89) 



(6.90) 



Note that the a s (m) corrections to the matching coefficients give rise to negligible next-to- next-to- 
leading order (NNLO) contributions in the ratios above. No further simplifications can be achieved 
at NLO without explicit assumptions on the counting. If the two states n and r are in the weak 
coupling regime, then 7£yw 3Sl ) = m(E n — E r ) |l(Jlj . In addition, the ratio of matrix elements in 
front of the rhs of l|6.89|l can be expressed in terms of the measured leptonic decay widths 



(Oi( 3 Si))„ _ T(V Q (nS) 



(Oi( 3 5i)) r T(V Q (rS)^e+e-) 



Img ee ( 3 6'i) E n - E r 
Im/ ee ( 3 5i) m 



(6.91) 



Img ee and Im/ ee are short distance matching coefficient which may be found in [3]. Ea. H6.9I|) 
and the expression for 7Z^, 3Si \ also hold if both n and r are in the strong coupling coupling- 
regime |S7| 11171 ITT] , but none of them does if one of the states is in the weak coupling regime 
and the other in the strong coupling regime. In the last case the NRQCD expression depends on 
five unknown parameters, which depend on n and r. If both n and r are in the strong coupling 
regime further simplifications occur. The matrix elements of the color octet NRQCD operators 
are proportional to the wave function at the origin times universal (bound state independent) non- 
perturbative parameters |57l 11171 ITT] (see appendix IU)l . Since (0i( 3 <Si))„ is also proportional to 
the wave function at the origin, the latter cancels in the ratios involved in l|6.90|l . Hence, TZ'q = 
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Figure 6.15: Background-subtracted CLEO data for the T(15) photon spectrum |48j . 



for the octet operators appearing in 
widths reads 



Then, the pNRQCD expression for the ratio of decay 



dz 
dT r 

dz 



T(V Q (nS) 



e T e" 



Imgee ( 3 Si) En - E r 
Im/ee ( 3 5i) m 



Cj [ 3 Si] (z) 1 
Ci [ 3 5i] (^) m 



(£„ - £ r ) (6.92) 



Therefore, in the strong coupling regime we can predict , using pNRQCD, the ratio of photon 
spectra at NLO (in the v 2 , (Aqcd/w) 2 |%71lll7| and a s (y/mAQCB) x v^QCd/ti m] expansion). 
On the other hand, if one of the states n is in the weak coupling regime, VJq will have a non-trivial 
dependence on the principal quantum number n and hence it is not expected to vanish. Therefore, 
expression l|6.92|) provides invaluable help for identifying the nature of heavy quarkonium states. If 
the two states are in the strong coupling regime, the ratio must follow the formula l|6.92|l : on the other 
hand, if (at least) one of the states is in the weak coupling regime the ratio is expected to deviate 
from H6.92fl . and should follow the general formula l|6.89|l . We illustrate the expected deviations 
in the plots (dashed curves) by assigning to the unknown IZs in l|6.89[l the value v 4 (v 2 ~ 0.1), 
according to the original NRQCD velocity scaling. 

We will use the recent data from CLEO |3H1 (which includes a very precise measurement of 
the Y(15) photon spectrum, as well as measurements of the Y(25) and T(35) photon spectra, see 
figures 16.151 16.161 and I6.17fl to check our predictions. In order to do the comparison we use the 
following procedure. First we efficiency correct the data (using the efficiencies modeled by CLEO). 
Then we perform the ratios 15/25, 15/35 and 25/35 (we add the errors of the different spectra in 
quadrature). Now we want to discern which of these ratios follow ea. H6.92|l and which ones deviate 
from it; to do that we fit ea. H6.92(l to each of the ratios leaving only the overall normalization as a free 
parameter (the experimental normalization is unknown). The fits are done in the central region, 
that is z £ [0.4,0.7], where ea- H6.92[) holds. A good (bad) x 2 obtained from the fit will indicate 
that the ratio does (not) follow the shape dictated by ea. (|6.92ll . In figures IB. 181 16.191 and 16. 2f01 we 
plot the ratios 15/25, 15/35 and 25/35 (respectively) together with ea. (|6.92|) and the estimate of 
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Figure 6.17: Background-subtracted CLEO data for the T(35) photon spectrum |48j . 
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Figure 6.18: Ratio of the T(1S) and T(2S) photon spectra. The points are obtained from the CLEO 
data 0Hj- The solid line is ea. (|6.92|) (overall normalization fitted), the dashed line is the estimate 
of (|6.89|) (see text). Agreement between the solid curve and the points in the central (unshaded) 
region would indicate that the two states are in the strong coupling regime. 

H6.89JI mentioned above (overall normalizations fitted for all curves, the number of d.o.f. is then 45). 
The figures show the spectra for z E [0.2, 1] for an easier visualization but remember that we are 
focusing in the central z region, denoted by the unshaded region in the plots. The theoretical errors 
due to higher orders in a s (m) and in the expansions below l|6.92f> are negligible with respect to the 
experimental ones. For the 1S/2S ratio we obtain a X 2 /d-0.f.|is/2S ~ 1-2, which corresponds to an 
18% CL. The errors for the T(3S) photon spectrum are considerably larger than those of the other 
two states (see figures Iq"Td1 16.161 and I6.17|) . this causes the ratios involving the 35* state to be less 
conclusive than the other one. In any case we obtain % 2 /d.o.f.| 1 5/ 3 5 ~ 0.9, which corresponds to 
a 68% CL, and x 2 /d-0.f.| 2 s/3s ~ 0.75, which corresponds to an 89% CL. Hence, the data disfavors 
T(15) in the strong coupling regime but is consistent with T(25*) and T(35) in it. 

In summary, using pNRQCD we have worked out a model-independent formula which involves 
the photon spectra of two heavy quarkonium states and holds at NLO in the strong coupling regime. 
When this formula is applied to the Upsilon system, current data indicate that the T(2S) and the 
T(3S) are consistent as states in the strong coupling regime 7 whereas the T(IS') in this regime is 
disfavor. A decrease of the current experimental errors for T(25*) and, specially, for the T(3S*) is 
necessary to confirm this indication. This is important, since it would validate the use of the formulas 
in |87l 11171 II lj . and others which may be derived in the future under the same assumptions, not 
only for the T(2S*) and T(35) but also for the Xb(2P )s, since their masses lie in between, as well as 
for their pseudoscalar partners. 



7 Y(25") also seems difficult to accomodate in a weak coupling picture in the analysis of the radiative transition 
T(2S) -+ Vbl \M- 




Figure 6.20: Same as fig. EH for T(25) and T(35) 



Chapter 7 

Conclusions / Overview 



In this thesis we have employed Effective Field Theory techniques to study the heavy quark sector 
of the Standard Model. We have focused in three different subjects. First, we have studied the 
singlet static QCD potential, employing potential Non-Relativistic QCD. With the help of that 
effective theory we have been able to determine the sub-leading infrared dependence of that static 
potential. Among other possible applications, this calculation will enter in the third order analysis 
of t — i production near threshold. An analysis which will be needed for a future e+ — e~ linear 
collider. After that we have studied an anomalous dimension in Soft-Collincar Effective Theory. 
That effective theory has very important applications in the field of £>-physics. A field which is of 
crucial importance for the indirect searches of new physics effects (through the study of CP violation 
and the CKM matrix). And finally we have studied the semi-inclusive radiative decays of heavy 
quarkonium to light particles, employing a combined use of potential Non-Relativistic QCD and 
Soft-Collinear Effective Theory. Viewed in retrospect, that process can be seen as a nice example 
on how a process is well described theoretically once one includes all the relevant degrees of freedom 
(in the effective theory) and uses a well defined power counting. When we have the radiative decay 
understood, it can be used to determine properties of the decaying heavy quarkonia, as we have also 
shown in the thesis. 



91 



Appendix A 

Definitions 



In this appendix we collect the definitions of some factors appearing throughout the thesis. 

"fE is the Eulcr constant -fE = 0.577216... . ((s) is the Riemann zeta function, with £(3) = 1.2021... 
. The Euler beta function is given by 

A.l Color factors 

The color factors for an SU (N c ) group are given by 

C f = ^ C A = N C T F = \ (A.2) 

A. 2 QCD beta function 

The strong coupling a s = g 2 /(Air) constant runs according to 

da s „ f „ a s . „ /a s \ 2 



The first coefficients of the QCD beta function are given by 

114 34 20 

A) = yCU - -T F n f ft = -C\ - -C A T F n f - AC f T F n f (A.4) 

where, here and throughout the thesis, rif is the number of light flavors. 
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Appendix B 

Feynman rules 

B.l pNRQCD 

The pNRQCD Lagrangian (|3.1U|) gives the position space rules for the vertices and propagators 
displayed in figure IB. II Feynman rules in ultrasoft momentum space are also useful. These are 
displayed in figure IB. 21 (additionally an insertion of a correction to the potential 6V in a singlet or 
octet propagator will give rise to a —iSV factor). 

B.2 SCET 

The Feynman rules that arise from the Lagrangian (|3.23|l are represented in figure IB.3I 
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Appendix B. Feynman rules 




Figure B.l: Propagators and vertices in pNRQCD in position space. We have displayed the rules at 
leading order in 1/m and order r in the multipole expansion. If one wants to perform a perturbative 
calculation these rules must be expanded in g. 



B.2. SCET 
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Figure B.2: Propagators and vertices in pNRQCD in ultrasoft momentum space. P M is the gluon 
incoming momentum. Dashed lines represent longitudinal gluons and springy lines transverse ones. 
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Figure B.3: Propagators and vertices arising from the SCET Lagrangian i|3.23|) . We have just 
displayed the interactions with one and two collinear gluons, although interactions with an arbitrary 
number of them are allowed. Dashed lines are collinear quarks, springy lines are ultrasoft gluons 
and springy lines with a line inside are collinear gluons. 



Appendix C 

NRQCD matrix elements in the 
strong coupling regime 



First we list the four-fcrmion NRQCD operators that appear in the subsequent formulas. 

Oi( 3 Si) 
OsCSo) 
8 ( 3 Si) 

V 1 ( 3 S 1 ) 



^crT a X ■ x f <xT> 

1 Vxx^-i^ 



2 

l r 

2 



D) 2 V + H.c. 
^•xM4d)V + H.c. 



Os( 3 Po) 



O s ( 3 P 2 ) 
Oem( 3 5i) 

PemCSo) 
~Pem( 3 Si) 



i) = 



^(iD)r x . x t(-iD)T^ 

i ^(-Id-o-^xxU-I D-<r)T> 

i ^(-f 5 xo-)r x • x f H D xo-)T a V 

V> f (-f D ^'^T^H 5 (VfljTV 
V' ctx| vac) (vac | X ct V 1 
1 " 
2 
1 
2 



■0^x1 vac) (vac j ( — -D 2 ) V-> + H.c. 



^<j-^|vac) ( vac lx <T 



D V + H.c. 



(C.l 
(C.2 
(C.3 

(C.4 
(C.5 
(C.6 
(C.7 

(C.8 

(C.9 
(CIO 

(C.ll 
(C.12 



In the strong coupling regime (Aqcd 3> E) the following factorized formulas can be derived for 
the NRQCD matrix elements. The analytic contributions in 1/m for some S- wave states (we just 
display here expressions involving S- wave states, since are the only ones really used in the thesis, 
see jSZ] for a complete list), up to corrections of (D(p 3 /m 3 x (Aq CD /m 2 , E/m)), are given by 



(V (n5)|Oi( 3 Si)|V (n5)> 1 / ,B = c j fl "°[° )|2 



E^ 2£ 3 



2£- 



(2:*) 



3m 2 3m 2 



(C.13) 
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Appendix C. NRQCD matrix elements in the strong coupling regime 



(V Q (nS)\0 EM e Sl )\V Q (nS)y/™ = C A ^>L ( 1 - ^ 2 -f + + ^ J (C.14) 

(V r Q (nS)|7'i( 3 5 1 )|Vg(n5)) 1/Tn = (Pg(n^|Pi( 1 So)|i > <j(nS)) 1/Tn = 
(VQ(n5)|P EM ( 3 5 1 )|F Q (r 1 5)) 1 / m - (P Q (n5)|P EM ( 1 5 )|P Q (r 1 5)) 1 /'" 

(VQ(n5)|0 8 ( a 5i)|Fg(n5)) 1 /« = (P Q (ri5)|0 8 ( 1 1 So)|PQ(r 1 5)) 1 /'" 

= Ci ^p(^f'j (CJ6) 

(V Q (nS)\O s eS )\V Q (nS))^ = (^WM^M^ 

o 

= c A ^?!f f (C A /2-C/)4BA 
27r \ 3m 2 / 

(^(n5)|O 8 ( 3 P J )|y o (nS0> 1/m _ <P Q (^)|0 8 ( 1 P 1 )|P Q (n5)) 1 /™ 



2 J + 1 9 

So 

2tt 



^^^(^■j (CJ8) 

(C.19) 



There are also non-analytic contributions in l/m. Up to corrections of order 0(p 3 /m 3 x Aqcd/iti x 
ma s / y/m Aqcd) they are given by (see f° r the complete list of known corrections) 

(V Q (nS)\0 1 ( 3 S 1 )\V Q (nS)} 1 /^ = (V Q (nS)\0 EM ( 3 S^VginS)) 1 /^ 

__ C ,K( 1 + 1 ^,^ (C , 0) 

In all those expressions R represents the radial part of the wave function, E the binding energy and 
all the £ and B are universal (bound state independent) non-perturbative parameters. 
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